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HAUPTAUFSÄTZE 


Über die Stabilität dreischichtiger Stäbe und Platten, deren 
mittlere aus einem Leichtstoff bestehende Schicht einen in 


Dickenrichtung veränderlichen Elastizitätsmodul hat. 1. 


Von F. Reutter in Karlsruhe. 


Bei Verbundstäben bzw. -platten mit zwei dünnen. Sperrholzplatten oder Blechen als Deckschichten 
und einer aus einem Leichtstoff bestehenden Füllschicht können zwei verschiedene Instabilitäten auftreten: 
das langwellige Bulerknicken des Gesamtstabes und ein kurzwelliges Kniltern, wobei jedes der beiden Bleche 
unabhängig vom andern beult. Diese Instabilitäten werden im folgenden untersucht unter der besonderen 
Annahme, daß die Füllschicht einen in Dickenrichtung veränderlichen E-Modul besitzt. 


The author discusses the properties of bars respectively plates composed of two thin plywood or metal 
sheets and a filling layer between them consisting of « material of small weight. With these bars, two different 
instabilities may occur; the Euler breaking of the whole bar with a long wawe-length and a crumpling with 
short wave-lengths. In the latter case, one sheet crumples independently from the other. T'he instabilities 
er lc on the presumption that the modul of elasticity of the filling layer is variable transverse to 
the. bar. 


Les barres ou les plaques composes, constitues par les plaques de bois ou de metal avec matiere legere de 
remplissage, peuvent donner lieu ü deux genres differents d’instabilite: la flexion transversale d’Euler ü grande 
longueur d’onde de la barre entiere et le froissement local de plaques superficielles independemment l’une de 
Vautre. Ces instabilites sont analyses par la suite, en admettant une variation du modul d’elastieite de la 

"matiere da remplissage dans le sens de l’epaisseur. 


B cocTaBHEIXx CTep’KHAX WIM INACTHHAX C ABYMA TOHKUMUH IOKPOBHEIMH CIIOAMU u3 
haHepsL mm Meranıa u CepmueBHHOÄ N3 JeTKorO MaTepnana MOTYT Habmparca Ba Pona 
HeyCTOÜyYHuBOTO PaBHOBecuA: HOMepeyHsHÄ Lporn6 Io Yünepy Bcero CTEPIKHA HIIH MECTHO® 
CMHHAHHE, IIPH KOTOPOM KaskıbIä N3 TIOKPOBHEIX CIIOeB BEITY4HBACTCTA He3saBuMO ONUH OT 
Aapyroro. Hu:ke päccMmaTpuBamTca 3Tu ABIEHHA HEYCTOÄYUUBOIO PaBHOBeCHH, IIPHHHMAA, 
YTO 3ANOAHAOINMH CHOA HMeeT B HaIpäaBlIeHHN TOJIHHBI TepeMeHHBÄ MOAyAIb YOPYTOCTH. 


Vor einiger Zeit ist der Vorschlag gemacht worden!), auf achsialen Druck beanspruchte 
Stäbe und quer- und längsgedrückte Platten im Leichtmetall- und Holzbau aus zwei verschiedenen 
Stoffen zu konstruieren. Zwei dünne Bleche oder Sperrholzplatten als Deckschichten sollten 
durch eine aus einem Schaumstoff bestehende leichte Füllschicht als Schubverband verbunden 
werden. Dieser Gedanke wird durch einen bekannten Sachverhalt aus der Festigkeitslehre nahe- 
gelegt: Die Festigkeit eines Stabes kann ohne Änderung der Querschnittsfläche durch geeignete 
Dimensionierung erhöht werden, indem man die Massen in möglichst großer Entfernung von den 
Hauptträgheitsachsen anordnet. Dies führte schon seit langem auf die bekannten Hohl- und 
Kastenquerschnitte. Allerdings ist der Verringerung der Wandstärke und Vergrößerung der 
lichten Weiten bei solchen Konstruktionen eine Grenze gesetzt, unterhalb der der Querschnitt 
dadurch instabil wird, daß an Stelle eines seitlichen Ausknickens des Stabes als Ganzes eine Wellung 
der Wände (Knittern) eintritt. 

; Wie die Hohlquerschnitte, so können auch Verbundstäbe zwei verschiedene Arten von 
Instabilitäten zeigen: Einmal kann der Stab als ganzes knicken (langwelliges Eulerknicken), 
zum anderen kann ein kurzwelliges Knittern eintreten, wobei jede der beiden Deckschichten 
unabhängig von der anderen beult. Welche von den beiden Knickmöglichkeiten sich einstellt, 
hängt von dem Verhältnis der Last zu den verschiedenen Bemessungsgrößen des Verbundstabes ab. 

In der vorliegenden Arbeit, die im Anschluß an eine Zusammenarbeit mit Herrn W. Flügge 
entstand, wird nun in Weiterentwicklung des obigen Gedankens das Verhalten von. Verbund- 
stäben und -platten untersucht, bei denen die Füllschicht nicht mehr homogen ist, sondern an den 
Rändern eine größere Dichte besitzt. Dann ist auch der Elastizitätsmodul der Füllschicht über 
die Dicke der Füllschicht hin in entsprechender Weise wie die Dichte veränderlich. Daher benötigt 
man für diese Untersuchungen die Elastizitätstheorie eines Mediums mit in Dickenrichtung ver- 
änderlichem E-Modul. 


ı) w.Flügge und K. Marguerre, Die optimale Knicklast eines Stabes, der aus zwei durch einen leichten 
Füllstoff verbundenen Blechen besteht. Luftf.-Forschg., Untersuchungen u. Mitteilungen (U. M.) 1360/1; K.Marguerre, 
Die optimale Beullast einer längsgedrückten gelenkig gelagerten Platte, die aus zwei durch einen leichten Füllstoff ver- 
bundenen Blechen besteht. Desgl. U.M. 1360/2. 


_ in der x- bzw. y-Richtung. Dann läßt sich mit m = 


[3 = 
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In Nr. 1 werden für eine elastische Halbebene und einen elastischen Parallelstreifen, dessen E 
E-Modul in der Dickenrichtung vom Rande ausgehend stark abnimmt, die Randverschiebungen 


unter der Annahme sinusförmiger Normalbelastung des Randes angegeben. Diese Normalbelastung 


wird durch den Druck der sinusförmig knitternden Bleche auf ihre elastische Bettung hervor- z 


gerufen. Die Füllschicht hat man sich als einen elastischen Parallelstreifen vorzustellen, dessen 


E-Modul beiderseits der Mittellinie von einem sehr kleinen Werte nach dem Rande zu stark ansteigt. _ 


Bei genügend kurzwelligem Knittern beulen die beiden. Deckschichten unabhängig voneinander 
so, daß jede von ihnen als auf einer elastischen Füllschicht von unendlicher Dicke (elastische Halb- 
ebene) gebettet angesehen werden kann. In Nr. 2 wird die „‚Bettungsziffer‘ bestimmt. Je kürzer 
die Wellenlänge, desto weniger unterscheiden sich die Werte der Bettungsziffer für Halbebene 
und Parallelstreifen. Die Füllschicht wird jetzt vorwiegend am Rande mechanisch beansprucht 
und die Deformationen dringen nur schwach in die Tiefe ein. In Nr.3 und 4 wird der auf Halb- 
ebene und Parallelstreifen elastisch gebettete achsial gedrückte Stab behandelt und seine kritische 


Stauchung bzw. Knitterlast ermittelt. In Nr.5 wird ein allgemeiner Weg zur Behandlung des 


Eulerknickens des Gesamtstabes angegeben. Die Differentialgleichung für die Verschiebungen 
erhält ein von der Knicklast abhängiges Störungsglied. Die Gleichgewichtsform des ausgeknickten 
Stabes läßt sich durch Iteration ermitteln, die Knicklast P selbst durch Nullsetzen einer Achter- 
determinante. Dieser sehr umständliche Weg kann vereinfacht werden, wenn man P näherungs- 
weise statt in die Gleichgewichts- in die Randbedingungen einführt. Weit einfacher ist freilich 
die Berechnung der Knicklast nach der Energiemethode, die im zweiten Teil dieser Abhandlung 
gegeben werden soll. 


1. Randverschiebungen für die elastische Halbebene und den Parallelstreifen 
unter sinusförmiger Normalbelastung bei exponentiell abnehmendem’ E-Modul. 


Wir behandeln zunächst kurz die Halbebene und den Parallelstreifen, deren E-Modul senk- 
recht zur Randlinie bzw. der beiderseitigen Berandung stark abnimmt. Die Randlinien seien sinus- 
förmiger Normalbelastung ausgesetzt. Es sei die Halbebene x > 0 
(Bild 1) gegeben, und wir machen für den E-Modul E(z) den 

Mey Ansatz 

ED) SS MEI ZEIT Free (1). 


Die Normalspannung am Rande x = 0 sei sinusförmig verteilt, 
und die Verschiebungskomponente parallel zur Randlinie soll ver- 
schwinden. Es seien v, v» die Komponenten der Verschiebung 
7 
i* 


Bild 1. 


‚21* Wellenlänge der sinusförmigen Normal- 
belastung, setzen: 
u(2,y) = p (x) sinmy +q (x) cos my 
v(%, y) = Ss (x) sin my +t(x) cos my h 


Man erhält p (x), q (x), s (x), &(x) als Integrale linearer homogener Differentialgleichungen mit 
konstanten Koeffizienten?). Unter den obigen Randbedingungen wird bei verschwindender Quer- 


1 
kontraktionszahl » = Fr (m Poissonsche Zahl): 


2.(2) = (a. FU) U ea (3), 
m 2. m m 
(x) = = +5 n) NL 2 et U NEE 
TR Te = 
ıa)=0, sa)=0 
mit 
ee 
g1=5 1% me NET TR SR (4a) 
und 
J 02 2 p 2 
201 EMRE ern 
do, 203 E, ’ a, — or 0? E, ee a a ee EN (4b). 


°) F.Reutter, Die periodisch belastete Halbebene und der Parallelstrei itin Di ic 
Elastizitätsmodul. Erscheint im Ing.-Arch. ton zoll In DIOKBNEIGHURNE VOEAUABEH OHREN 


— 
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Insbesondere ergibt sich dann für die maximale Randverschiebung 
| k fe 

yser een 2 

ET ae 


0 Hr nn aa og a ee 
?(0) ee : 
| 23 TR | 
Anschließend behandeln wireinen Parallelstreifen von der Breite 25 (Bild 2). Der Z-Modul gehorche 
dem Gesetz | 


(E„m : E-Modul in der Mittellinie x = 0), es findet also ein Sprung von in der Mittellinie statt. 


Beide Ränder seien sinusförmiger Normalbelastung unterworfen. Wir unterscheiden den asymme- 


X 
Bild 2. 2 Bild 3a (asymm.). Bild 3b (symm.). 
trischen und den symmetrischen Belastungs- und Verschiebungsfall (Bild 3a bzw. 3b). Die Rand- 
punkte sollen keine Horizontalverschiebung erleiden. Dann lauten die Randbedingungen: 
o,(+b, y) = P,sinmy . ..:. (7a), %(—b,y)=F P,sinmy, . .. ..(7b) 
ee (8a), DE a RR . (8b). 


Das obere Vorzeichen in Gl. (7b) gilt für den symmetrischen, das untere für den asymmetrischen 
Fall. 4 
Unter diesen Annahmen kann man setzen?) 


u,y)=p(e)-sinmy, v(,y)=tla)cosmy .!. 2... (9) 
p(x) und t(x) genügen linearen Dgln. mit konstanten Koeffizienten und zwar für 
a) 2>0:. 2” +29" + (k2— 2m?) p” — 2m: kp +m?p = \ Nr (10). 
b) x<0: p"’ —2Kkp"" + (k?— 2m?) p’ +2m?kp Hmtp—0 
Ihre Integrale sind mit 
Be EIERN ee 
=> er tm, = = y# ee (11), 


a) Pu (2) = (a +a,2)e&® +(m, Faza)en | 


RE (12), 
b) Bal) = (An + Ay ame HlAs + Asa)ene | 
2 m m m 
a = u )« [0 au| ee + 
) [74 ) Mm 0? 1 01 0 01 ; 
2 m 
+ (+ as +” as +" a,alenr, 
m 05 02 
9 hr PR EP ee ee. 3 (13) 
p) =|(- )A A " Ayalemer 4 
) tu(®) alas ln, n 
2 m m m 
+4 +2)4. 42442 4,2]eos. 
BT, 01 01 
Zur Festlegung der insgesamt 8 Integrationskonstanten dienen die 4 Gln. (7), (8) in der Form: 
‚ ES Pr y. Fr ir, * 
RO: a re N ne RE U) 
1.0) =.0,, DENE DT (8*). 
Hierzu kommen in der Streifenmitte die Bedingungen 
Bo Bo el | 3 a 


P:(0)=Pr41(0), »% (0)=p4 (0). J 
®2) F. Reutter,l.c« 
1* 


| Reutter, 3.8 d 


= We +ya,, 
R a, a +4na, 
ri 20 se Der ST 
RR Me +22). 
NER a ya 
ei mtb m? DR 
ES: y= eteon (DiE uprT‘ + er NR 
Be ET ran (Mi = a) 
Be = — end Zi’ e=e Pr] Er 
mer med DE ) 
Br Re n= eare)b ve Ver mim: ’ N >. 
e wo das obere Vorzeichen für den symmetrischen, das untere für den symmetrischen de Fi Br 
E . so ergeben sich für den asymmetrischen Fall die Gleichungen PN | Ei 3 
E: : WER \ 0,= 4A, a,= — 4, u„=A, u =—4;- ne un. u 
A @,, @, genügen dem System 2 


mit 


m 
» 


Be PB: 
aß +0), +m— a7 +1), = aa —6 Ar ee 
Bere—etß+e)a +Bein +30: —el Y)a= gj%— 30° Be: 
Für den symmetrischen Fall erhält man dagegen 
= 34, 9, = 4, = Än AS a en) 
und a, = —a,. 


So lassen sich a,, @,, a; aus den drei Gleichungen: a | = 

. E a +Y4s Rh 

? q SET OT | 
er] EEE N - (20) Fe 
1=Öd—e aa+NQz . x 

ea. —201, +2, =0. Bi 

bestimmen. . 


Dann sind im symmetrischen wie im ee Fall p (x) und t(x) und somit auch 
u(z, y), v(z, y) und hieraus die Spannungen O4 Oys Ta, bestimmbar. 


2. Die Bettungsziffer für Halbebene und Parallelstreifen. 
. Die sinusförmige Normalbelastung 
Fo=PR my 
bedingt beim elastischen Halbraum die Randverschiebung 
(0, y) = p (0) sin my (21), 
bei der elastischen Halbebene dagegen 
u(b, y) = p (b) sin my (22). 
Normalbelastung und Randverschiebung koppeln wir gemäß dem bekannten von Zimmer- 
mann formulierten Ansatz) durch die Relation 


Py)=—aulD, ya ee ne » (28) 
bzw. 
Py)=-—x' u(b, Y) Er U ee ie (24). 
x heißt die ‚„‚Bettungsziffer“, Dieser Ansatz erweist sich als brauchbar für die kurzwelligen Vor- 
gänge (Knittern), während sich außerhalb dieses Bereiches auch hier seine in der Literatur mehr- 
fach diskutierte Problematik zeigt®). 
‚ Nach Gl. (4b) wird jetzt 


4 


‘) Zimmermann, Die Berechnung des Eisenbahnoberbaues. Berlin 1888;u.Hayaschi, Theorie des Trägers 
auf elastischer Unterlage. Berlin 1921. 
) K.Wieghardt, Überden Balken aufnachgiebiger Unterlage. Z. angew. Math. Mech. Bd.2 (1922), S. 165—184; 


K. Marguerre, Über den Träger auf elast. Unterlage. Z. angew. Math. Mech. Bd. 17 (1937), S. 224—231. 


TB 


> B,= - B-Modul der Füllschicht am Rande, E= B-Modul der Deckschi i 


R a L x _m*, A (Pr) sind dimensionslos. 
wird - 


\ 
4 


2 0 

5 2 Al ee Fer - 
(B%) ist ein von ee Biockdicke t und dem Modulverhältnis u unab hi ängiges Maß für die 5 
Bettungsziffer. Die Abhängigkeit der Bettungsziffer x bzw. f} (82) von ß, 1* und k* ist aus Tafel1 
Bichtlich. 2 Ziel der Untersuehungen dieser Nr. ist ein Vergleich der Beitunesallet, für die } nr 


Tatel 1.7482) rei: 


< „213| 0,003689| 0,01 


0,0496 1 0,2222 | 0,2981 | 0,4714 | 0,6667 | 0,8164 0,9428 1,3333 


20.3 0,00884 | 0,0268 0,1632 | 0,2462 | 0,4185 | 0,6126 | 0,7641 

2.0,3- = 0,00463 | 0,08889 | 0,1626 | 0,3231 | 0,5204 | 0,6630 | 0,7863 | 0,9973 | "— 
0,5 10,000202| 0,00133 _ | 0,04887 | 0,1070 |] 0,2501 | 0,4504 | 0,5720 | 0,6933 | 0,9000 | 1,1007. 
0,6307 | — 0.00072 | 0.03363 | 0,0811 \ 0,2092 | 0,3799 | 0.5161 | 0,6365 | 0.8411 | 1.0148 

10  10,000027 | 0,00019 | 0,01317 | 0,0386 | 0,1282 | 0,2676 | 0,3921 | 0,5000 | 0,6916 | 0,8579 
2,00 en 0,00005.  |.0,00222 | 0,00797 | 0,0373 | 0,1051 | 0,1804 | 0,2555 | 0,4000 | 0,5353 
30 | — .| 0,00001 | 0,00069 | 0,00266 | 0,0141 | 0,0469 | 0,0885 | 0,1346 | 0,2330 | 0,3333 
5,0 — | 0,0000017 | 0,00016 | 0,00061 | 0,00385 | 0,0133 | 0,0277 | 0,0453 | 0,0885 0,1386 


“ Halbebene mit der für den Parallelstreifen. Dazu müssen wir eine zu f, (2?) analoge Funktion 
für den Parallelstreifen ermitteln. Nach Gl. (22) und (24) wird 


SR ) r Hl t 

N. n (9) = B,uM3  peb)-wis | 

wo » (b) gemäß den Gln. (11), (18), (20) einzusetzen ist. Die Bettungsziffer ist jetzt von dem Ver- 
hältnis 


. (27), 


us —bhr 


FE _ abhängig. 3 #R 

er Gemäß den GIn.(11), (15), (16) ist 2(b) proportional zu d; x und f, (ß?) sind dagegen umgekehrt 
_ proportional zu b. | 
- Die Tafeln 2 und 3 enthalten die Werte von: f* ( (89 = ; (ß°) für den symmetrischen 


bzw. asymmetrischen Fall in Abhängigkeit von bk* und 52 ß2, ae Tafeli den Grenzfall 
-b> ©o wiedergibt. 


be pr 
231 0,014656 | 0,04 0,4 0,8 2,0 4,0 6,0 8,0 12,0 


Tafel 2: f#(ß?) (Parallelstreifen, symmetrischer Fall). 


0 0,60861 | 0,60876 0,61352 0,62551 0,67782 | 0,78030 | 0,88815 | 0,98989 | 1,17753 
0,2 0,60362 | 0,59084 0,59649 0,60422 0,64698 | 0,74255 | 0,84249 | 0,94061 | 1,12176 
1,0 . | 0,52687 | 0,52570 0,51366 0,50423 0,52120 | 0,58920 | 0,67017 | 0,76620 | 0,93174 
2,0 0,39398 | 0,39276 0,38210 0,37541 0,34342 | 0,42566 | 0,49230 | 0,56884 | 0,74669 
4,0 = == 0,17236 0,16998 0,17299 | 0,20341 | 0,25047 | 0,30663 | 0,42781 
6,0 — Fi 0,064479 0,06352 0,06835 | 0,08919 | 0,12079 | 0,15987 | 0,24805 
0,0 — 0,00625 0,005982 0,00651 0,00792 | 0,01792 | 0,03134 | 0,04859 | 0,09027 


0.8884 | 1.1001 | 120785. 


Pe a Air eh en TBSEN (4 ch Du Ya fe REIT. ah 
h . F « Y. u 
x 


> u SR) 


J 4, i a RUER 2 u 2 
j > + L Ye R EmnT F = A ok 1] eT Je 
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Tafel 3: f (ß) (Parallelstreifen, asymmetrischer Fl. 
0  [0,003033 | 0,00se01 | 0,07861 | 0,15084 | 0,34056 | 0,57055 | 0,75174 | 0,89768 | 
0,2 [0.000162 | 0.002999 | 0,04665 | 0,10962 | 0,28569 | 0,51688 | 0,69708 | 0,84247 | 
1.0 |0.000028 , 0.000192 | 0,01842 | 0,04517 | 0,14683 | 0,34143 | 0,54363 354 | 
2.0 [0.000011 | 0.000073 | 0,00605 | 0,02079 | 0,08930 | 0,22372 | 0,35294 | 
4,0 Ex = 0.00170 | 0,00635 | 0,03206 | 0,09604 | 0,16943 | 
6,0 — Fe 0,00079 | 0,00246 | 0,01366 | '0,04511 | 0,08582 | 0,13208 | 
10,0 — | 0,000038 | 0,00034 | 0,00058 | 0,00253 | 0,01321 | 0,02735 | 0,04522 


Die Diskussion der numerischen Ergebnisse (mansetze am besten inden Tafeln2und3b=2cm) 
zeigt: Je größer ß?, um so geringer werden die Unterschiede zwischen den Zahlenwe der 
Tafeln 1, 2 und 3. Für große 2, d.h. kleine Wellenlängen kann also die Halbebene an Stelle 
des Parallelstreifens gesetzt werden, ohne daß ein nennenswerter Fehler begangen wird. Sehr stark 
unterscheiden sich auch für kleine ß? die Bettungsziffern im symmetrischen und asymmetrischen 
Fall beim Parallelstreifen; die asymmetrischen Werte liegen viel tiefer. Wir kommen darauf noch 
bei der numerischen Behandlung des Eulerknickens im 2. Teil dieser Arbeit zurück. 


3. Die Knitterlast für ein auf der elastischen Halbebene gebettetes Blech. \ 

Wir betrachten jetzt ein einzelnes Blech, das mit dem Rande der Füllschicht verbunden und ER 
daher auf dieser elastisch gebettet ist (Bild 4). Die Dgl. für die Durchbringung % = % (y) des ela- F 
stisch gebetteten Stabes, der an beiden Enden y= (0 und y=1 a 
dem achsialen Druck P unterworfen ist, lautet: . . 


EJ,wW"" +PW’ +ua=0 3 a (28) 
Bild 4. Dabei bedeutet x wieder die „‚Bettungsziffer“. 
3 
Mit dem Ausdruck (25) für x und J, = 5 lautet jetzt die Dgl. (28): m 
zu 122, 24 07 E, ar ie . 
U — Et U E,(m® +20?) % uU= 0 N er (29). 


Wir fordern.die Randbedingungen der gelenkigen Lagerung und behandeln die Dgl. (29), 
indem wir % nach den Eigenfunktionen des achsial gedrückten, gelenkig gelagerten (nicht elastisch 
gebetteten) Stabes von der Länge / entwickeln: 


| %(y) = 2 ön SIn.my Ye rt Dr Se (30). 
Mit 
SAN EX 0 E23 240? E, 
My = ] = Ir und B = — E, (m? +2 0°) ln = (30a) 
erhält man dann für P: 
E, £ 4 4 
= IDPRIER (7 4 Bl ) Re Se Be ee (31). 


Für diesen Wert von P ist eine von % = 0 verschiedene Gleichgewichtsform des Bleches möglich. 
Inc 4 7 A = R I Stablänge 
Diese ‚‚Knitterlast‘‘ P hängt also nur vom Verhältnis !’* = — — SE mia | =W ä “ni 
' ’ g | nis N Wellenzahl Wellenlänge, nicht 
von der Stablänge absolut ab. Will man auf die Stablänge beziehen, so muß man eine diskrete 


Kurven schar, d.h. für jede Wellenzahl N eine andere Kurve für die Knitterlast angeben. 


Für k=0 wird, =—m 
und 
E\®n 22,1* 
limP= u tr « E 
ser jan TS ee Re Se (32); 
hieraus folgt als Grenzwert der kritischen Stauchung € — P/Eit: 
5 ti? n® 2.5, 1* 
lim. 2= im Fe er 
k->0 k>02,t 12179 3. Hunt. ©. 27 Besen (33). 


P ist die Knitterlast für ein Blech; die Knitterlast für den Stab, d.h. beid i 
Einsetzen des Wertes von B (Gl. 30a): ee 


(20: +2,) (319). 


°) Siehe z.B. Biezeno-Grammel, Techn. Dynamik 8,521, 


Tr 
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Als Maß für die Knitterwellenlänge benutzen wir wieder die in Nr. 2 eingeführte Größe ß = ae 
dann kann die kritische Stauchung ® für ein Blech in der Form 


DE 2tuo? 
el, RR ER N A) 
20, P° Yoe 
geschrieben werden oder mit Hilfe des Ausdrucks (25) für die Bettungsziffer auch 
een 34 
e= D E,ß® EL Der re ( a) 


Man zeigt leicht, daß & proportional zu „2/3 ist. 


Die sehr kleine Wellenlänge /* wird sich so einstellen, daß &ein Minimum wird. Wir schreiben 
€ mittels 


= (0) — koı) : 
als Funktion von og, und setzen zur Bestimmung des Minimums a2 —(, 


d 
Es wird & 
ee t2 kt2 2tuo 
= 2 Er 1 
e=- Bin a7z Bot ate ea (35) 
und hat einen Extremwert, wenn e, eine Wurzel der Gl. 5. Grades: 
19 34 4u 19 5 ka: A 
Bir Ringer +28) etz ra — (75 +36 ®) DE else 


Für k=0 wird o? +20, also imo = y—4n. 
k—0 


Für k £0 behandeln wir die Gleichung am einfachsten 
mittels der Substitution o, = z%k. Dann ergeben sich die Null- 
stellen als Schnittpunkte einer festen Kurve 5. Grades mit 
einer Schar von Parabeln 2. Grades mit gemeinsamer Achse 
und festen Nullstellen (Bild 5): 


19 34 19 5 1 4 u 
VIER | Er ET NIERTT a Er Bu 2 
2 6° + a In +9 2 18 3m 322) (36a). 


Es zeigt sich, daß unabhängig vom k-Wert eine posi- 


Er 
tive einfache Nullstelle o, in der Umgebung von Be 


vorhan- 


/0,25081 u 
Weitere reelle Nullstellen sind nurfür k < IE 


denundzwar sind diese negativ. 
/0,25081 u 
t 


Bild 5. 


Für k,= 


strichelte Parabel berührt eben die Kurve 5. Ordnung. Für k < k„sind je zwei einfache negative 


liegt eine doppelte Nullstelle mit 2 = — 0,9624 vor, d.h. die ge- 


kt 
Nullstellen, für k > k, keine reellen negativen Nullstellen vorhanden G = m) > 


BER oe k? RICH 
Nun ist nach Gl. (4a) o, eine negative Zahl. Zudem bedingt, gleichgültig ob Y vs + Pr positiv 
oder negativ, 0, — n- einen imaginären Wert von ß. Es ergibt sich also: Nimmt der Z-Modul der 


Füllschicht (unendliche Halbebene) nach dem Ansatz E (x) = E, ekzah, so gibt es ein Minimum 
Ne en w 


Jedoch ist die Berechnung der Knitterlast unter Zugrundelegung der unendlichen Halb- 
ebene für k*-Werte, die wesentlich von Null verschieden sind, nicht ohne weiteres zulässig. Wir 
ermitteln die Wellenlänge, die zum Minimum der Knitterlast gehört, für k = 0 und für k, für die 
Annahmen t = 0,067 cm, u = 1,668 - 10-* und d=2cm und erhalten 


Km24t, Emtt<2b, 


der kritischen Knitterstauchung (mit reeller Wellenlänge) nur fü 


er. Bi RT 


“sehr bedeutend sind. 
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aber doch in der Größenordnung 2b. Die Länge der Wellen des knitternden 
im letzteren Falle schon mit der Dicke der Füllstoffschicht vergleichbar. In 


auch die tiefer in der Füllschicht liegenden Teilchen bei der een 3 


spielen. Daher muß zur Behandlung des Knittervorganges die Rechnung für den E 
Parallelstreifen mit ausgeführt werden. Allerdings zeigt die Durchführung, daß die tatsächlichen 
Abweichungen zwischen Parallelstreifen und Halbebene im Gebiet kleiner Wellenlängen n 


- 


4. Die Knitterlast für ein auf dem Parallelstreifen elastisch gebettetes 


Die geschlossene Angabe der kritischen Stauchung & für den Parallelstreifen, dessen E-Modul 
beiderseits der Mittellinie exponentiell ansteigt, würde auf einen sehr verwickelten Ausdruck 
führen. Wir berechnen daher &nach Gl]. (34a) mit Hilfe der in Nr. 2 schon angegebenen numerischen 
Werte der Bettungsziffer und können dazu die Knitterlast auf Grund der Beziehung 


Bü me 
Pe m tr 


= 
bestimmen. Nach Gl. (34 a) ist & wiederum proportional zu „*'®. Die Abhängigkeit der kritischen 
Stauchung up von ß? und k* bei einem Stab von der Dicke 2b = 4 cm für den symmetrischen 


Fall ©, den asymmetrischen Fall ® und die Vergleichswerte für den Parallelstreifen ® zeigt 
Tafel 4. Für andere b-Werte lassen sich die Tafeln leicht aus den Tafeln 2 und 3 berechnen;eine 
dimensionslose Tafel wie für die Bettungsziffern läßt sich indessen hier nicht angeben. Doch 

ist mit der Kenntnis der &-Werte für ein endliches 5b und den Grenzfall b— » (Halbebene) das 
qualitative Verhalten von & für verschiedene b-Werte ausreichend gekennzeichnet. Wie in Nr. 2 

für die Bettungsziffer zeigt sich auch hier: Je größer ß?, um so geringer werden die Unterschiede 

der kritischen Stauchung für Halbebene und Parallelstreifen. Die Resultate der Halbebene stimmen 

bis auf eine geringe Senkung der Knitterlast weitgehend mit denen für den Parallelstreifen (asymme- 
trischer Fall) überein (Tafel 4, Zeile ®). 


2 
- 
a 


IE ©® Parallelstreifen symmetrisch. 
Tafel 4: —575 für | ® Parallelstreifen asymmetrisch. 
u 


® Halbebene. b=2cm. 
I(em)| 63,00 38,24 12,09 8,56 5,41 3,82 3,12 2,70 2,21 
I* 2 
nn 0,003689 0,01 0,1 0,2 0,5 1,0 1,5 2,0 3,0 
u | 
® 164,98 60,88 6,144 3,144 1,397 0,8636 0,7171 0,6616 0,6425 
0 ® 0,8205 0,8610 | 0,7945 0,7708 0,7228 0,6539 0,6262 0,6155 0,6274 
® 10,98 4,960 2,231 1,507 0,9845 0,7500 0,6693 0,6381 0,6349 
® 163,63 59,09 5,973 3,038 1,336 0,8259 0,6867 0,6370 0,6239 
0,1 ® 0,04431 | 0,3008 0,4749 0,5650 0,6131 0,6002 0,5897 0,5879 0,6074 
® 2,397 2,680 1,612 1,247 0,8777 0,6959 0,6343 0,6109 0,6167 
PIE, NENERESPERENN EOBERRRE ENT SUR, USE SCH SE DEREN RETTEN EN ©. . 
® 142,82 52,57 | 5,145 2,538 1,084 0,6725 0,5718 0,5498 | 0,5606 
0,5 ® 0,00802 | 0,02005 | 0,1926 0,2425 0,3353 0,4248 0,4874 0,4949 0,5423 
® 0,05526 | 0,1338 | 0,4971 | 0,5514 | 0,5417: | 0,5358 | 0,5063 | 0,5133 | 0.5500 
® 106,80 39,28 3,829 1,894 0,7285 0,5090 0,4532 0,4511 0,4989 
1,0 ® 0,00322 | 0,00816 | 0,06046 | 0,1206 0,2203 0,3071 0,3603 0,4008 0,4792 
® 0,00755 | 0,01993 | 0,1400 | 0,2098 | 0,2977 | 0,3509 | 0,3863 | 0,4167 | 0.4805 
® ame —_ 1,732 0,8666 0,3877 0,2867 0,2920 0,3200 0,3926 
2,0 ® = = 0,02533 | 0,04844 | 0,1058 0,1794 0,2380 0,2895 0,3808 
® = == 0,03052 | 0,05651 | 0,1163 0,1884 0,2453 0,2944 0,3833 
® ur — | 0,6531 | 0,3343 | 0,1784 | 0,1725 | 0,2055 | 0,2466 | 0,3327 
3,0 ® — 0,01626 | 0,02896 | 0,06901 | 0,1284 0,1822 0,2327 0,3269 
® — = 0,01518 | 0,02999 | 0,06987 | 0,1302 0,1840 0,2340 0,3277 
® 0,6260 | 0,06816 | 0,03422 | 0,05752 | 0,1013 0,1459 0,1910 0,2801 
5,0 ® == 0,00469 | 0,01176 | 0,01956 | 0,04675 | 0,09654 | 0,1432 0,1893 0,2798 
® 


0,00100 | 0,00991 | 0,01974 | 0,04937 | 0,09663 | 0,1435 0,1893 0,2796 
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; 5. Die Knieklast des Verbundstabes. 

Zur Behandlung des Knickproblems kann der Stab nicht als schlank angesehen werden, 
doch könnte die Knicklast bei konstantem E-Modul der Füllschicht mittels der durch die Schub- 
korrektur modifizierten Euler-Formel berechnet werden (s. II. Teil). Für die Füllschicht 
von veränderlichem E-Modul schlagen wir hier zunächst ein strenges Verfahren ein. 

Wir betrachten jetzt den gesamten Verbundstab aus zwei dün- 
nen Deckschichten und der Füllschicht; er sie gelenkig gelagert und 
in Richtung seiner Achse einer mittels einer Backe-an greifenden 
Druckkraft P unterworfen. In Bild 6 ist ein Element des unter der 
Einwirkung von P ausgeknickten Stabes dargestellt. Sind o, 0, 
0, die Komponenten des Spannungstensors, so ist o, nicht mehr 
klein von erster Ordnung, es kann aber in der Form 


%,=0o0W)-+e, Bild 6. 


angesetzt werden, wo o\ konstant und o, von 1. Ordnung klein ist. Für die Verschiebungs- 
komponenten gilt jetzt der Ansatz (2) und mit 


E 
Hr Eu: 0, Bi). 7,> en (4, +9,), eo) = : ‚ P Knickkraft, 
erhält man für die Funktion ? (x) in Gl. (2) im Falle x > O0 die Dgl. 
[2220 270 123 fi m? Pe" [2 
p”’ +2kp" +(® — 2m?) p' —2m?kp Be — kp! —2m?p). (37). 
m 


Für P=0 ist sie mit der Dgl. (10) Nr. 1 identisch. Wir betrachten daher die rechte Seite von 
Gl. (37) als Störungsglied und versuchen die Gleichung durch Iteration zu lösen, indem wir das 


Integral» (Gl. 12a) von Gl. (10a) als nullte Näherung benützen und in die rechte Seite von . 
Gl. (37) einsetzen. 


Mit 3% k? 3k k2 
zei — 2 ee er 2 
a a tat”: d E+yotn 
wird als erste Näherung für ein Partikulärintegral der inhomogenen Gl. (37) angesetzt: 
m (br +Hbra)e "Hs Hr a)e"”. 


m 


Mit bF = = b, wird dann 9, = pr) + pı eine 1. Näherung für die allgemeine Lösung der Gl. (37): 


m 


P 2 = 
2,= la ee} Fe a N ee 
£3 +7 rar ko ) 


PR ® x 
a a nen TE a he 
+le +, rear at )| 


m 


Nach Ausführung des Koeffizientenvergleichs erhält man: 


Dee rmfn te. (a,Q, +a,Q, a, 0:0 ar eh) 
mit 
(gs merke, we Ener? 16. 
. (39 
EFF EL en P 20, (k?— m?) —m?k Be 
RZ 16. bKög, a 16 bK?o, 
Ganz entsprechend ergibt sich 
ee a +8 +98,2) +e’(, 8, +, +8,80). . ... (40), 
wobei 
| —m P (eı + k) (ea — 2%) PN? —0:(0ı +k)lee—2k) 
u ER 01 \0ı 2 5 2. 10ı 
ee un. | = ln) 64m bE R2 


Amer 16 mbk? 
® 


a Bar nn Hartıhı,. |? \V ale he ram 
a ta \4mb 16 m b k? Ber 64 m b2 R? 


rt DR E R 
ee TED B- 


NEE 
ER 


dreischichtiger 
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Be Pp LT GERDRDS 
ee? Tomb 


FESTEN 


PP a2 BO +MREm) 
| m 2. P ,—1 8&—2h(@+2h , 2 —M(2—m) IR 
Al ers ReLzEl ET TE > 16 m b le 
Bi. pP (ala Fb at A K% 
Be: | ar) ge 64 m b? k? 0, a 


Analog ergeben sich q (x) und s (x). Die Funktionen p, t, s und q können allerdings nicht 
für den ganzen Stab einheitlich angesetzt werden, sondern es gelten in den einzelnen Bereiche 
die aus Bild 7 ersichtlichen Werte, wobei 


PL (x) == (& +0, x) ASeh + (&5 +0; x) Pas 


i,(@) = (a, + e) PERS +( +2 0. +22) 


Pıy (X) = (a, +01 a) Marie +(0) +0) a) Pi 


’ 3 ’ ’ =: 2 3 ’ , mE 
(a) = (a Me o)e x +(«; 6% +0, 2) e 5 


während für ?r, fr die Werte 9,, Ta (Gl. 38, 40), für Pyır, 

£;rr die analogen Funktionen ?,, t4 zu setzen sind. 
Entsprechend bezeichnen wir auch die @,, 0, Tzy | 

mit o} usw. (Der E-Modul der Deckschicht wird kon- 

stant und so groß angenommen, daß die rechte Seite der 

Dgl. (37) vernachlässigt werden kann, also die Lösung der 

homogenen Gl. ausreicht.) Wir suchen die Knicklast P 

unter folgenden Randbedingungen: 


Oberer Rand der Deckschicht («=b +b): 


SE ee See er (41a) rn, 6e+9=0... u rue (41b), 
Unterer Rand der Deckschicht: 

RAIL IT U) Se (41e) &.(6) et (B) 2.5 2 200 ar De (41d) 

N ee ie (4le) 22 Heid): - +»... 0.. (41f). 


Indem man die Faktoren der sin- und cos-Glieder je für sich vergleicht, zerfallen die GIn. (41) 
von vornherein in 2 völlig identischen Systeme, von denen das eine nur p und £, das andere nur 
qund senthält. Es genügt also, mit den Funktionen p und tallein weiterzuarbeiten. Für die Mitte 
der Füllstoffschicht (x = 0) gelten die Bedingungen 


ET (A1g) OB OT RR (41h). 


(Das Knicken entspricht immer dem asymmetrischen Fall — Nr.2 —.) Hier ist in der Streifen- 
mitte o,—=0 und v=0, was mit (41g und h) identisch ist, während 


Pr (0) = Py7(0) und pi, (0) = Pin (0) 
sich automatisch aus der Bedingung der Asymmetrie ergibt. \ 
Hierzu kommen die entsprechenden 6 Bedingungen für den unteren Blech- und: Füllstoff- 

rand, die nach Einsetzen der Funktionen Pyyp +++, Zy die Größen A,, a enthalten. Gemäß 
Gl. (17) ist in diesen 6 Gln. und in (41a—h) zu setzen: 

a,= 4, a4,=— 4, 

ER ,=—o 
Dadurch werden die 6 Gln. für den unteren Rand mit den Gln..(41 a) bis (41 f) identisch 


und zusammen mit (41 g) und (41 h) verbleibt ein System von 8 linearen homogenen Gin. 
für die 8 Unbekannten a,, @, Ay, Ay, &y &, &, &. 


0 


In den Bedingungen (41 a, b) und (4le,f) ist in den I-Größen die Querkontraktionszahl 
der Deckschicht mitzuberücksichtigen, also 


«+? 
JE 


vy der E, 


U 
%,=#&, Ay ng, +%,) 


dabei ist: 


% 2 mt H-[mi(u-i) +] 
Denn 2m 


-69%+3 


| Se Dr | 1 
een: 
EN) BEREITET 
: all #9)... - ern, 


Praktisch ist mehr als eine Iteration nicht-ausführbar und auch diese nur mit sehr erheb- 


lichem Rechenaufwand. Für Z= 63 cm ergeben sich aber mit einer Iteration noch keine brauch- 
baren Resultate. Dies ist verständlich, da eine zahlenmäßige Überprüfung für die rechte Seite 
der Dgl. (37) dann Werte in der Größenordnung 1 ergibt, für 1= 5,41 cm dagegen in der Größen- 
ordnung 102. In diesem Falle erhält man für k* = 0,5 und u = 16,68.10°, E, = 7: 10° kg em? 
und » = 0,3 den Wert P= 177,75kg, der recht gut zu dem gemäß der Klementaren Energie- 


. methode nach Timoshenko erhaltenen Wert (s. II. Teil) paßt. 


. Eine wesentliche Vereinfachung ergibt sich, wenn man von vornherein nicht von En 
Dgl. (37) ausgeht, sondern von Gl. (10), die Gleichgewichtsbedingungen also nicht am verformten 
Stabelement aufstellt, sondern die Verformungsgrößen auf die Randbedingungen wirft. Unter 


{ 


Beibehaltung von pr, &i, ?ıv, tiv verwenden wir 9% Pa ta ta — Gin. (12) und (13) — für pır, - 


Prir, bır, kru- 
Es ist dann die Normalspannung am oberen Rand der Deckschicht gleich der durch den 


Achsialdruck erzeugten Querbelastung q zu setzen. Dies führt im Falle » = 0 auf die Gln. — wir 
lassen künftig das identische System in qg und s wiederum weg —: 
Bo) —mPpb+)=0 B-abt+y)—mPgb+)=0.. . (4a). 
Zweitens soll die Schubspannung am oberen Rande verschwinden: 
RIEF Ele (42h) 


(und dieselbe Gl. in qg und s, die wir künftig nicht mehr angeben). 

Am unteren Rand, wo Deckschicht und Füllschicht zusammenstoßen, verlangen wir Über- 
einstimmung von Vertikal- und Horizontalverschiebungen, Normal- und Schubspannungen, d.h. 
u) um. , ad, Bo) Mob)... eb 

E, (mpı (db) +1 (6)) = E, (mp (db) Hi (Ö)) 200. (421), 
während in der Mitte x—=0 wieder die GIn. (41g) und (41h) gelten. Jetzt geht P nur linear in die 


WR ee” are 


"Knickdeterminante ein, und.diese läßt sich so. 


BP a 8 
PB r 


Äh 


=D 
=; > 
2 
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ı - Pr 5 a 5 e = ee 
minanten erscheint. Wir geben diese in Gl. (43) an, indem wir gleichze tig ( 
der Deckschicht mit berücksichtigen: En. 5 


1 Ham? mt (et mt 41) (+2 eure) % 
A ee in REM: Be 
Seen —y er er; 
SH, In 9.0 0a Ion | 
item t| 2 mn ee 
ar N) +0 (& _& euren) 
I+» 02 0 
ne (l-—e-atob+ : 
g amd mt (em — 1) 2 Alde.li+ 2? -aren) 
| tz (ei nn) | 
| ae De Rlm(1+22e-mtene) + zumlı a +00) 2 
al De er user ae Te 02(, ,0 BER: 
Ei EHE) a 
mi + 92 e(a+ oo) mb (1 -& e-(uı+e)b Mr 
01 1 E 
en (545) — (+ ateo) 
2 
E+2 (0002) 2(1 8 tat) 
(em +1) | 3 (em 1) — 2 mt BEE re 
7 03 (1 —eato)dL 
(e2mt +1) \ f \ Abe. (1 En N ea+ oo) 
tzei-e)| 
aim(ı rn ara) möm(1 _& aa) 
a a a = ' ; ie 
| I (5 7 5 e) (e1=e3)| 55 (1 u a) 


Es ist: u=u(l +»). 

Im Falle k*—=0 vereinfacht sie sich erheblich. Setzt man an Stelle der Bedingung 
Tyy(b +t)=0 die (in Nr.1 für v (b) eingeführte) Bedingung v (b +1) = 0, so ist die Gl. (42b) 
abzuändern und jetzt schreibt sich P sogar als Quotient zweier Dreierdeterminanten. Es wurden 
hierzu folgende Zahlenrechnungen ausgeführt: 


i=63 cm k*=0, v=0, „= 16,68° 105: P = 196,258 


I= 855m H=0, ve, P = 224,97 kg 
*=0, v=03 P= 247,22 kg 
*=05 v=03 P = 154,95 kg, 


während sich mit der Bedingung v (b +1) = 0 mit » = 0,3 sogar ergab: 
P(k*=0): P (k* = 0,5) = 356,84: 255,20. 


Die Ergebnisse zeigen, daß die im II..Teil mit der Energiemethode gewonnenen Knicklasten 
auch im Falle des exponentiell veränderlichen Z-Moduls eher zu niedrig als zu hoch gegriffen sind. 
Hierauf einige Stichproben zu machen, war neben der Darlegung der allgemeinen mathematischen 
Entwicklung der Zweck dieser Nummer. 


Eingegangen: 24. Juni 47. 


der laminaren Grenzschichten. | 
Teil II. 


Von = er in Freiburg i. Br. 


m 

glich, a, aha daher eine we ee der’ rg auf einige wenige Zahlen 
rte für o nicht notwendig macht, näher untersucht werden. Dazu A zuvor ein Wort Apr den Hi 
a raktischen Wert dieser Rechnung gesagt werden. _ Ri dh 


= Die Theorie ist wegen der Annahme U,= U,(y) von Karese auf die Vera 
0 <1/2% (siehe (5.4)) beschränkt. Zum anderen muß die i in Zusammenhang mit (5.19) gestellte 
allgemeine Forderung, daß nämlich die Dicke der Zone, in welcher die äußere Stromlinienwellig- 
keit außerhalb der Grenzschicht abklingt, groß ist gegenüber der Dicke der Grenzschicht‘ selbst, 
_ erfüllt sein, d.h. es muß o Re? (L)> 1 sein. Da Re(L) praktisch jedenfalls auf Werte unter 106 
beschränkt bleibt, damit die ungestörte Strömung laminar verlaufen kann, muß alles in allem 
105 < <o<1/2zsein, so daß o = 4/2 L auf Werte in der Nachbarschaft des kleinen Zahlenwerts 
10 unseres Beispiels aus $5 beschränkt bleibt. Das wird also auch für die folgende Theorie 
gelten, die demnach gegenüber dem Beispiel des $ 5 zunächst nur insoweit einen Fortschritt bringen = 
kann, als sie die analytische Abhängigkeit der Ergebnisse von o in dem verhältnismäßig engen 
zulässigen o-Intervall um a — 107? aufzeigt. Solange man jedoch keine umfassendere Theorie e- 
- sitzt, wird man die folgenden Ergebnisse wenigstens größenordnungsmäßig auch extrapolatorisch 
für Werte von o und Re verwenden, die den obigen Voraussetzungen nicht mehr entsprechen. 
"Es erscheint daher berechtigt, die folgenden Ausführungen zu bringen, die einen tieferen Einblick 
In die analytischen Zusammenhänge und größenordnungsmäßig brauchbare Regeln bieten. Sie 
wollen aber wie das Vorangehende nur als erster Vorstoß über die Berechnung numerischer Einzel- 
- fälle hinaus gewertet werden. 


Wir gehen davon aus, daß das Blasiussche Plattengrenzschichtprofil (5.6) 


2 4 2 
U en en E = - 


(& = 0,33206) in Wandnähe wegen = in y= . verschwindenden Krümmung durch das Hncade 
Profil 5 
/ 


EEE ri 


gut oröximiert werden kann. Noch für y=1 (also bis zum”Wandabstand ö,, siehe Bild 1) 

liegt die Abweichung unter 1% des strengen Wertes. In hinreichender Wandnähe kann man daher - 

für (4.3) genähert 
; . - ayf, -afıtegy =1 


RES a0, anal = 


schreiben, oder zusammengefaßt für die komplexe Funktion 
Be | ky)=hly) +in(y) mit i= 1) 
- oh ayh re MN. ee (6.3). 


Das allgemeine Integral der linearen Differentialgleichung (6.3) lautet: 


mtr Bm l+oy 2.22... 


wo Zaj3 die Algemeipe Zylinderfunktion der Ordnung 2/3 bedeutet: 
ZEHN ine... un (6.42) 
(©, Os, © Integrationskonstanten, Be ’Hankelsche Funktionen 1.bzw.2. Art). Schreibt man 
2= Yolo y layer e (6.5) 


Be REN ES 
Y en, ne, = : . m, TR ie ne %, Gr hi 
in N A AS | Ding a : k Bet, Kate, 
- (reelle Wurzel), so wird A als Funktion von 2: 2 


Na je) + aa (re) + 


Nimmt man vorübergehend an, daß die Näherung U, = a y nicht nur in Wandr 
jn0 Sy <1), sondern in 0 Sy < ® gilt, und verlangt man neben h (0) = #’ (0) =! 

_  h für y—- endlich bleibt, so werden — siehe unten — durch diese von o unabhäng en 
bedingungen die Integrationskonstanten c,, csund cz in (6.4) bereits eindeutig festgelegt. Nto)=0 
Jäßt sich wegen U,— © für y—- © nicht mehr erreichen.) Die Störungsamplitude A (alp ı 


De 


und g,) wird also eine universelle Funktion von z. Es gibt eine Zahl z= z,, die ein Maß für die 
Dicke ö, der „‚Störungsgrenzschicht‘‘ abgibt derart, daß die Amplitude der Störungsgeschwindig- 
keit, also A’ (z), in 2 = 2, praktisch (mit vorgebbarer Genauigkeit) ihren konstanten asymptotischen — 
5 _ Werterreicht hat. Es ist wegen z = z, = Zahlsomit y= , = Zahl+ o\/?, also, daymit ö,dimen- 
 sionslos gemacht ist, die Länge ö, —= Zahl - 0/3 ö,. Wegen o = A/2z.L ist somit 


\ 


P2 
4 I 


% YA f | ne Air 
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0 8, YvL/UO, ist ein Maß für die Dicke der wirklichen ungestörten Plattengrenzschicht. Istnur 


_ A/L hinreichend klein, so ist also die Dicke ö, der Störungsgrenzschicht beliebig klein gegenüber be. 
der Dicke der ungestörten Grenzschicht, insbesondere so klein, daß die Zone von der Wand bis 

Br zum Wandabstand 6, innerhalb jener Zone liegt, in der in praktisch ausreichender Näherung wirk- 
000 lich O,=ay gesetzt werden kann. Für y>1 muß die Krümmung des Profils U,(y) berück- 
4 sichtigt werden. Hier schließt sich die mit ö, von o abhängige jeweilige asymptotische Lösung 
an. Die wesentlichen relativ großen Störungsbeträge liegen innerhalb der Zone bis zum Wand- 
abstand ö,. Hier allein (bis h’ (z) seinen konstanten asymptotischen Wert erreicht, siehe Glieder 
mit o in Gl. (6.2)) wirken sich die Reibungskräfte maßgeblich auf den Störungsverlauf aus. (Die 
Bo Störung greift natürlich, der vorgegebenen äußeren Geschwindigkeit U,+ ecos zentsprechend, 
Be“ bis in die äußere Potentialströmung hinein.) Als „‚Störungsgrenzschicht‘‘ bezeichnen wir hier so 
jene wandnahe Zone, in der durch die Störungen hervorgerufene Reibungskräfte wirksam werden. 
Hier erfolgt die relativ stärkste Störung, während außerhalb der Störungsgrenzschicht die Rei- 
bung für den Verlauf der Störungen nur noch gemäß dem Verlauf von U, (y) schwach ins Gewicht 
fällt. Für Einzelheiten sei auf die folgende Rechnung verwiesen. 


Als wichtiges erstes Ergebnis dieser Betrachtung ist somit auszusagen: Ist das Verhält- 


r nis der Störungswellenlänge A zur überströmten Wandlänge Z nur hinreichend 
: klein, so wird die „Störungsgrenzschicht‘“, in der die Reibungskräfte sich auf = 
en :. | den Störungsverlauf wesentlich auswirken, nur eine relativ sehr dünne wand- 
nahe Zone („Untergrenzschicht‘) innerhalb der ungestörten Plattengrenzschicht a 
BD einnehmen. Ihre Dicke wird proportional ValE klein im Verhältnis zur unge- 
Ber störten Grenzschichtdicke. as 


; Für Wandnähe wurden oben (siehe (4.4)) Reihenentwicklungen nach Potenzen von y ange- 

„ setzt. Unter der Annahme U,= &y wird die Ermittlung der Beiwerte dieser Reihen besonders 
einfach. Überdies gestattet die — für das praktische Rechnen unmittelbar wenig geeignete — 
Lösung (6.4), die in den Reihenentwicklungen zunächst freibleibenden Parameter ß, und y; leicht 
zu. berechnen. Dies sollim folgenden geschehen. Es werde also nunmehr weiterhin der Grenzfall 
sehr kleiner o-Werte angenommen. 


: Um aus dem allgemeinen Integral (6.4) das hier gesuchte partikuläre Integral auszusondern, 
müssen einmal an der Wand die Randbedingungen A (0) = h’(0) = 0 erfüllt werden, zum anderen 
muß bei Zugrundelegung des Profils U,(y)=ayinO0 <y< ® für wachsende y zumindest ver- 


langt werden, daß h’(y) endlich bleibt. Setzen wir der Einfachheit halber die untere Integral- 
grenze y, in (6.4) zu 


fest, so ergibt schon die letzte Forderung in (6.4) mit (6.4a) 


N  .. (6.7),9 


wie man ohne weiteres aus dem bekannten asymptotischen Verhalten der Hankelschen Funk- 
tionen folgert. Es ist dabei angenommen, daß o so klein ist, daß die asymptotischen Werte in der 


in welcher in SWirklichkeik ge ER 4 et ten En Bakisch erreicht worden. Es | = x 


ird sich zeigen (siehe unten), daß etwa o < 10? sein muß. 


Er; Die Bestimmung von C, erfolgt allein auf Grund der Forderung k(0) = 0. Es ist nach leichter 

6: arme unter ne der Ba Nücking von Ay : 
me: ;: ORZESCEE ö EINER 

E.* Fe on man ‚Ve N /2 |) _—t— g1/3 32/3 6 all sin2/3), 

3 WE 2 WR 0 7) S, 

Brise | | 

E- ER G=—oniB nr, )a- e2”il3sin Bu) RN = ER 6.8). | Km 


E: bekannt sind, 


(0) = ar.) e2ils [343 ro) a a = (6.9) 
somit wegen ky)—=f(y)+ ig,(9) | 
| ; Ba 2 () ER 
Be. “= j j ; FERIEN 2 3 a 
E- ee zur (?) Br N Re Ba 
& s Für das numerische Beispiel o= 10 des $5 ergab sich (siehe (5.8)) fa = 21,1, 9, = —37,3. 


gleich ß, = 21,3, y; = — 37,0. Während in $5 die Werte von ß, und y, mit Hilfe des Anschlusses 
. der wandnahen Lösungen an die asymptotischen Lösungen mit einem zulässigen Rechenfehler von 
etwa 1% gewonnen wurden, ergeben sie sich hier aus den als Funktionen von z = y!?(a/o) 
universellen wandnahen Lösungen unter der Voraussetzung hinreichend kleiner o-Werte. Man 
erkennt, daß für o = 10 die Übereinstimmung bereits erzielt ist. Wir werden weiter unten 
erkennen, daß in der Tat o = 10°? bereits einigermaßen in den Bereich der hier zulässigen kleinen 
o-Werte fällt, so daß eine gute Näherung erwartet werden mußte, womit sich zugleich eine gegen- 
seitige Bestätigung der verschieden gearteten Rechnungen von $5 und $6 ergibt. 


= Man kann vermöge (6.10) bereits den Verlauf der Wandschubspannung angeben. Es ist ja 
in den gewählten Dimensionslosen (siehe (5.14)) 


zn (2,0)=&+ e{P,cosx + yssin «} 
E- 9 Im“. 0 ne (6.11) 
E RE” al1-+e —Lleos (+ 3) 


‚Nach dieser Theorie der 1. Ordnung ar eist also der Reibungswiderstand der Wandwellenlänge A 
der gewellten Wand gleich jenem der gleichen Länge der glatten Wand. Die Wandschubspannung 
schwankt periodisch um den Wert der Wandschubspannung der ebenen Platte (im betrachteten 
Intervall einiger Wellenlängen, in welchem unsere approximativen Ansätze gelten). Erst die hier 
nicht berücksichtigten quadratischen Glieder in e bringen eine Abweichung der Reibungswider- 
stände der ebenen und der gewellten Platte. Diese Abweichung wird also in zweiter Ordnung 
bezüglich € klein. Wir kommen am Schluß dieses Paragraphen auf die Frage des Widerstandes 
— zum Reibungswiderstand kommt noch der Druckwiderstand hinzu — ausführlicher zurück. 
In erster Ordnung bezüglich &, also für hinreichend kleine Welligkeiten, gilt nach (6.11): Die 
Wandschubbspannung läuft der Schwankung des äußeren Druckes pP — m, = 
— 0 D2ecos x» mit einer Phasenverschiebung von 2/3 = 120° nach. (In dem numerischen 
Beispiel o—= 102 des $5 ergab sich 119,50 als Phasenverschiebung.) Der äußeren ‚wirksamen‘ 
Welle der Stromlinien am Rande der Grenzschicht (siehe (5.17)) läuft die Wandschubspannung 
entsprechend in der Phase um x/3 = 60° voraus. 

Nach (6.9) ist «&/|h’’ (0)| = 0,16764 02°. Diese Zahl gibt nach (6.11) denjenigen e-Wert an, 
für den bei von O0 an en eerstmals 0U/dy an der Wand verschwindet (und zwar in 
x —= 120°, also 120° hinter dem Minimum, 30° hinter dem Wendepunkt und 60° vor dem Maximum 


E:.. - » Ebenfalls mit Hilfe der Reihenentwicklung von H%h, ergibt sich nunmehr, da C, und c 


Nach (6.10) wird für o—= 10”? wegen & = 0,33206, I'(1/3) = 2,6789, I'(5/3) = 0,90275 zum Ver- - | 
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des äußeren Drucks). Als Ablösungskriterium ergibt sich somit die Aussage: Ist 2 
| 5> 0,1080. zu) ST Ra . (6.12), Br 
der Grenzschicht ein. Dabei erfolgt Rückströmung n 


so tritt Ablösung 
das ecos («+ 3) < —0,168 0? ist, also innerhalb eines 


Wandnähe in jenem x-Intervall, für 
zu x —= 120° symmetrischen Intervalls. Leider gibt e erst ein Maß für die Amplitude der wirk- 


samen Welle am äußeren Ran 5 > 
ermittelt werden (siehe unten). Das numerische Beispiel o = 10? des $5 liefert zum V 


wegen &xrit = 0,00775 in guter Übereinstimmung mit (6.12) als kritische Zahl 2/o®® = 0,167, was 
wiederum als Kontrolle der dortigen Rechnungen dienen kann, da der kleine Wert o—= 10, 
wie schon im voraus bemerkt wurde, den Voraussetzungen für die Gültigkeit der vorliegenden 
Theorie einigermaßen entspricht (siehe unten). 

Das zahlenmäßige Resultat des einfachen Gesetzes (6.12) beansprucht natürlich nur größen- 
ordnungsmäßige Richtigkeit. Man bedenke, daß ‚Störungen‘, die zu einer Ablösung führen, im 
Rahmen einer linearisierten Störungstheorie kaum in mehr als größenordnungsmäßiger Näherung 
erfaßt werden dürften. Gerade hier bieten nun die von Quick und Schröder [9] numerisch 
berechneten Beispiele eine willkommene Kontrollmöglichkeit. Von den vier dort mitgeteilten 
Beispielen scheiden die beiden letzten für diesen Zweck sicher aus, da sie Fälle betreffen, die auf 


“o-Werte von ungefähr 0,3 bzw. 0,4 führen, während sich die vorliegende Theorie, siehe (5.4), auf 
co <1/2r beschränkt. In den dortigen Beispielen wird zunächst als Anlauf ein ebener Plattenteil. 


der Tiefe L, überströmt, an das sich von hier an eine gewellte Platte anschließt. Im ersten Beispiel 
ist © = 0,005 und A/L, = 0,072, im zweiten e = 0,01, A/L, = 0,72 (im dritten Beispiel e = 0,025 

A]L, = 0,72, im vierten Beispiel e= 0,05, A/L,= 0,72). Die Länge Z unserer obigen Theorie 
gibt definitionsgemäß eine mittlere Länge an in dem Intervall einiger Wellenlängen, in dem die 
Strömung durch die gemachten Ansätze approximiert wird. Betrachten wir die Umgebung jener 
Wandstellen, wo inden Quick-Schröderschen Beispielen Ablösung zuerst eintritt, so ist 
sinngemäß L als überströmte Länge bis zur Ablösungsstelle zu wählen. Im ersten der genannten 
Beispiele entnimmt man den Quick-Schröderschen Abbildungen Z/Z, = 1,26, im zweiten 
Beispiel 4,30 (dagegen schon Ablösung bei 1,40 bzw. 1,31 im dritten bzw. vierten Beispiel). Es wird 
also im ersten Beispiel 2 Z/A= ao! = 110, somit eo?!®? = 0,115. Das ist bereits eine größen- 
ordnungsmäßig erfreuliche Übereinstimmung mit unserer Theorie, die den größeren Zahlenwert 
0,168 für hinreichend kleines ao liefert. Das zweite Quick-Schrödersche Beispiel ergibt 
o! = 37,5, also einen schon etwas kleinen Wert für den Vergleich mit unserer Theorie. Erstaun- 
lich ist aber, daß als kritischer Zahlenwert eo’? hieraus praktisch unverändert gegenüber dem 
obigen Wert 0,115 nunmehr 0,112 resultiert! Es dürfte demnach unsere Zahl 0,168 der Theorie 

die ja eigentlich nur für sehr kleine e und o anwendbar ist, etwas zu groß sein das allgemeine 
Ablösungsgesetz & - 02? = Zahl scheint jedoch weitgehend brauchbar zu sein. (Soweit ich gegen- 
wärtig voraussehen kann, bringt eine Theorie der 2. Ordnung in eeine Verkleinerung der Zahl 0 168 

aber sie dürfte noch immer merklich größer als 0,11 bleiben.) tr 

Um Aussagen über die Verdrängungsdicke der Grenzschicht und it ü 

der Wand zu erhalten, die zu der angesetzten periodischen Störung der ee 
führt, muß nun erst der Gesamtverlauf von f, und g, ermittelt werden. Dies soll zunächst in ade 
wandnahen Zone durchgeführt werden, wo U, durch das lineare Profil x y hinreichend a eis 

miert wird. In dieser Zone sind f, und g, universelle Funktionen der Veränderlichen z — 1 (a 
Be A von fi und 9, nach Potenzen von z, die aus (6.4b) aluelehee Be: 
ayanen te entstehen, zeigen hier besonders einfache Koeffizientengesetze. Setzt man 


Be Bi a nn 6.13 
1 @W)= —PrAz(e) + Yadı (2) Ira B;(2) \ 
mit den universellen Konstanten ß,,’, definiert durch 
2 2ni 
B. +iy Yu [07 AB: R 1 3. R 5 
2 Ya= = (B,+ %y,) — —2r(3)e 343 & re) er ee (6.14), 


also 


B, = 2,0655, NY SRG 20 vi (6.14a), 


de der Grenzschicht. Die wirkliche Wand muß wie in $5 erstnoch 


m 222 5 Dugk Su 265 


zu DR 
te 7: I 13° u 


RE 1E N: 141 201° 


; si 7 Ant 7. 10: 13: m —4:7- 10- 13-16- mt SR 
#5 Are —2.5:8- mtr 8: al 5:8- 11: 14:17:20: Bee 
Er — De +2. 5-8 11: 5-8-11- 14: 17: + 


Bei der numerischen Auswertung wurden die BR Ah Glieder berücksichtigt. Sie reichen RE 
aus, um nach (6.13) f, (2), 9, (z) und deren Ableitungen nach z von z = 0 an bis zu jenen z-Werten 
zu berechnen, für die fund g, praktisch fortan linear von z abhängen, df,/dz und dg,/dz also kon. 
‚stante Werte (ihre asymptotischen Werte bei Zugrundelegung des Profils U,=ayin0Z yZ oo) 
annehmen. Dies ist beietwa 2=3 erreicht, wie man aus dem in Bild 11 dargestellten Verlauf‘. na 


vn 


en ee dr 


Br 


% 


N DR 


Bild 11. Verlauf der Störungsfunktionen f, (2), g,(z) und ihrer Ableitungen nach z. 


= es Funktionen und ihrer Ableitungen erkennt. Man ermittelt aus diesem Bilde, daß für den asym- 
ee erreichten linearen Verlauf 


fı= 3,4 2 —3,32 
= — 2,00 2 + 0,27 ] 


ist. (Die Koeffizienten von fin (6.16) sind mit einem Fehler von maximal 1/,%/,, die von g mit 
maximal 2°/, ermittelt.) 

| - - Nungilt die Approximation der ungestörten Grenzschicht U‘, (y) durch ein lineares Profil & y 
in guter Näherung für,y <1 (siehe oben), also ist unsere Theorie der obigen wandnahen Lösungen 
für 2 < (a/o)V/® anwendbar. Andererseits muß z bis etwa zuWerten > 3 wachsen können, damit 
das durch (6.16) gekennzeichnete asymptotische Verhalten der Störungsfunktionen und ihrer ersten 
Ableitungen nach 2 in dieser wandnahen Zone mit 1°/, Abweichung erreicht wird. (Diese Annahme, 
daß die asymptotischen Werte bereits in dieser Schicht erreicht werden, liegt schon dem Ergebnis 
(6.7) zugrunde.) Aus diesen beiden Forderungen folgt mit « = 0,33 größenordnungsmäßig, daß 


co <10%2 


sein muß. Es bestätigt sich also, daß unser Beispiel aus $5 gerade noch von der Theorie des $6 
erfaßt wird, worauf schon die verschiedenen zahlenmäßigen Vergleiche hingewiesen haben. (Für 


2 
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grobe Näherungen könnte man sich schon mit z < 1,5 bis 2 zufrieden geben (siehe Bild 11), dart 


E erwarten.) be - 
ür o-Werte, die 10-2 etwas überschreiten, noch grobe Näherungsaussagen 2.) > 
ns ER —L d.h. z > (x/o)Y®, besitzen wir die asymptotische Lösung (4.9) 


A=4 U,+h—eBJ(D,), 5 
= BU, +0oAJ(U)+oJ(h) 


Die von o freien Glieder geben die 


Per ee" . 4 . i (4.9). 


‚reibungslose‘ Außenlösung, die 

“] Glieder mit o erfassen den Rei- 
ii bungseinfluß in erster Näherung für 

kleine o (siehe $4). Daß der Rei- 

et * pungseinfluß außerhalb der „Stö- 
9001 rungsgrenzschicht‘‘ nicht mehr we- 

9002 sentlich ins Gewicht fallen soll, be- 

= ©» deutet — und dies wird auch die 


folgende Rechnung bestätigen —, 
daß die letztgenannten o-Glieder 
—n 5 577 relativ klein bleiben. 

Die in (4.9) rechts auftreten- 
den Funktionen wurden für „,„=1 
’ bei der Auswertung des Beispiels 
oo  a—=10*von $5 bereits numerisch 
«os ermittelt (siehe Bild 2). A und B 

berechnen sich als Funktionen von 
2 co vermöge 


"50 


Bild 12. f,(y) und g,(y) für verschiedene Werte von o (wobei die u ( f U ) 
kleineren Werte nur theoretisches Interesse besitzen). — WuıYToly-l 


—= 10,43 (a/o)Y®— 10,07, 
B= (g/U )y=1ı 
—= —6,06 (a/o)Y3 + 0,8 


6.17). 


(Siehe hierzu (6.16) und setze . 
z= (a/o)\3sowie U ,[y= 1] =0,3298 
ein. Dabei ist wieder (a/o)/? 23 
vorausgesetzt, siehe oben.) In Bild 
12 und 13sind f, und g, bzw. fi und 
gı für verschiedene o-Werte als 
Funktionen von % aufgetragen. 
(Striche bedeuten Ableitung nach 
y.) Im Hinblick auf die einleiten- 
den Bemerkungen dieses Paragra- 
phen über die zulässige Größen- 
ordnung vono haben die kleineren 
o-Werte nur theoretisches Interesse. 
Durch (6.17) ist der stetige An- 
schluß von f, und g, in ya =1 ge- 
währleistet. fi und gi machen, wie 
Bild13. %(y) und 9 (y) für verschiedene Werte von o. die numerische Rechnung zeigt, 
(Die kleineren o-Werte beanspruchen nur theoretisches Interesse, in y=1 einen kleinen Sprung, 
siehe Text.) der allerdings durchweg so klein 
ist, daß er in der Auftragung 
von Bild 13 kaum merklich auffallen würde. Er ist von der Größenordnung der numerischen 
Ungenauigkeit in den Koeffizienten von (6.16). Die Kurven wurden hier entsprechend geringfügig 
graphisch geglättet. Für y >1 bringen die Reibungsglieder in der oben angeführten Außen- 
lösung (4.9) keine nennenswerten Beiträge. Dort, wo die Kurven des Bildes 13 für die Ge- 
schwindigkeitsamplituden f, und g; der Störung mehr oder weniger scharf umbiegen, hört der 
Reibungseinfluß auf den Störungsverlauf praktisch auf. (Trägt man übrigens zur Kontrolle die 


Kurven für «= 10”?in Bild 3 ein, so ergibt sich wieder gute Übereinstimmung mit den Rech- 
nungen des $5.) 
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Für y— © erhält man nach (5.11) wegen (6.17) 


[3 . 
im (A—y)= —6,4-+0,170+ he ® {10,43 —1,27 0), 


VER 


; 1/3 ; 
lim g, = 0,824 0,450 + (2) {6,06 — 2,18 0}, 


yon 


‚und 'also wegen o < 10°? (also (&/o)!® = 3,2): ‚ohne wesentlichen Genauigkeitsverlust weniger 
schwerfällig: 


ER PART . 13 
lim (h —y) = —6,24 + 10,43 04 (2) 
lim g, = 0,82 — 6,06 (=) 


(Die vernachlässigten Glieder gehen gerade aus den o-Gliedern in (4.9) hervor, die den Reibungs- 
einfluß auf die Außenlösung erfassen.) 


Die Verdrängungsdicke ö* der gestörten Grenzschicht ergibt sich (siehe. (5.20)) zu 
ö* = 65 —eÖ,D(o) cos (2 + %) 


mit öo = 1,72 ö, (Verdrängungsdicke der ungestörten Grenzschicht) 

und De)=[®+b]%, wy=—bla ..... BE HB LIN 
& 1/3 

= a= 4,52 — 10,43 (2) 


1/3 


b= —0,82 + 6,06 3 


Die Phasenkonstante y erweist sich für o < 10”? als nur schwach von o abhängig. (Von Re(L) 
ist sie unabhängig.) Für o = 10 ist y = 32,5°, und y nimmt mit abnehmendem o monoton ab, 
“ umim theoretischen Grenzfall o—0 gegen y = 30,2° zu konvergieren. Praktisch wird es genügen, 
mit » = 30° zu rechnen. 


Analog zu (5.25) erhält man als Wandstromlinie : 


Yin (a) = TE cos ae en, 
0 k 


wobei 

Ay= £ö,[(o Re” —a)? + 52]? 

ee 
(Re= Re(L)—= U,L]/v). Die Länge Ay, die Amplitude der Wandwelle, ist neben der Wellen- 


länge A als gegeben anzusehen. Daraus berechnet sich dann die Amplitude A,= e/ß der wirk- 
samen Welle am äußeren Rande der Grenzschicht nach (6.20a) vermöge 


4A, o Reli 
| TB ee ae (6.21) 
Schreibt man 6* = 6$ — Ag cos (© + y), so ist nach (6.19) und (6.20 a) ? 
Sr 22 He 1/2 
As Er re (6.22) 


Aw  |(o Rei? —a)2 + b 


Im Bild 14 und 15 sind die Ergebnisse (6.21) und (6.22) aufgetragen, wobei allerdings wegen der 
einleitend in diesem Paragraphen betonten Einschränkungen weite Teile der Kurvenbilder nur 
theoretisches und bestenfalls größenordnungsmäßig praktisches Interesse haben. Sie sollen in 
erster Linie die allgemeine Tendenz bei variablem A und Z deutlich machen und sprechen für sich. 


I* 


a Ce ern & ya 


“A > de 
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Bild 16 gibt die zugehörigen Werte der Phasenverschiebung p an. Sie ist für kleine o, wie. las YS 
Bild zeigt, praktisch wieder gleich = 30°. (Da in diesen Bildern Re und o mit Z gebildet.werden, % 
deutet man wachsendes o — A/2r L zweckmäßig als wachsende Wellenlänge bei festem Ort L.) R 


—_—— [) 
ooot 0002 be JR 71 Be 7, EEE |. ET, ao 002 Go a ae 
Bild 14. Verhältnis der Amplituden von wirk- Bild 15. Verhältnis der Amplituden von Ver- 
samer und wirklicher Welle. . drängungsdickenschwankung und Wandwelle. 


Das Ablösungskriterium (6.12) ergibt Ablösung, wenn die kritische Amplitude der Wand- 
wellung ; 


A yunıe = 168 0218 d, [lo Re? — a? +2... .... > 
. oder wenn die maximale Steigung $ Ay der Wandkontur den kritischen Wert 
(B Aw) 9168 07"? Re-V2[(o Re? —a® +BJR ....... (6.23) 


überschreitet. (Es ist 86, = 1/o Rel!2.) Es sei daran erinnert, daß eine genauere Theorie statt 
des Koeffizienten 0,168 eine etwas kleinere Zahl liefern wird (siehe oben). In Bild 17 ist (BP Aw)yrit | 
dargestellt. Es gilt bezüglich der Gültigkeit auch hier die zu Bild 14 bis 16 gemachte Bemerkung. | 
Insbesondere sind ’hier die Zahlenwerte durchweg nur größenordnungsmäßig zu verstehen in dem- 

selben Maße, wie die Zahl 0,168 unsicher ist. 


o MT —— — 
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Bild 16. Phasenverschiebung zwischen der Bild 17. Kritische Maximalsteigung der Wandkontur 


Schwankung der äußeren Geschwindigkeit 
Kr ee bei deren Überschreitung Ablösung erfolgt. 


Zum Schluß noch ein Wort über den Widerstand der welligen Wand. Es wurde schon 
oben betont, daß der Reibungswiderstand pro Wellenlänge sich erst in der 2. Ordnung bezüglich & 
gegenüber dem entsprechenden Wert bei der ebenen Platte ändert. Dieser Unterschied konnte 
daher in der bisher dargestellten Theorie der 1. Ordnung in & nicht erfaßt werden. In der allge- 
meinen Fourier-Reihe für die Stromfunktion 


Pe N Re BSACHO) cos 9x + b,(y) sin v e) 


wurden bisher nur die in e linearen Glieder a, = ef, (y) und , = eg, ( y) berücksichtigt. Das für 
die obengenannte Anderung des Reibungswiderstandes einer Wellenlänge A maßgebende Glied 
nullter Ordnung a,/2 erweist sich als in e quadratisch. Die Rechnung, die an späterer Stelle bei 
der Darstellung der Theorie 2. Ordnung wiedergegeben werden soll, liefert das Ergebnis, daß der 


Reibungswiderstand W,(4; Z) einer Wellenlänge A der welligen Wand in der Umge- 


bung der Stelle Z kleiner ist als der Reibungswiderstand W7, (4; L) des entsprechenden Teil- 
stücks einer ebenen Platte, und zwar wird 


Wr (A; 12)= Wr (4:2). 1 ae ee (6.24), 


1 (6.19) definiert ist 
ors ist z.B. im Falle von Bi 
e etwa 0,93. Es mache 


de 


 widerstand noch ein Druckwiderstand hinzu. Auch der Druckwiderstand ist von der 
EI 2. Ordnung ing. Er soll im folgenden berechnet werden. Zunächst soll jedoch noch der Reibungs- . 
FAN en der einseitig überströmten ebenen Platte der Länge L und der Breite B angegeben 
werden [1]: Fr | d8 SE 


: ER DT Wr(D)=e7'BL'4a: Re"W2(]) ar ee a (052 
A ee j > ? E > n DT RAR j. 
.. Er DerBeitrag desPlatienstücks von Z bis L+ 4 zum Reibungswiderstand ist also wegen Wr (Z)"L!2 
MED WEI OD Wed ApL] ee 
wer = Weimar] 
= Somit ist nach: (6.24) = | Be; 
E- 5 ü y R be& 
= - Wr; D)=MW;(L | ee en RE TEEN 
Br Wa (1) = Walze 1 - a | Ba 
E::. - Der Druckwiderstand der welligen Wand berechnet sich aus dem vorgegebenen Druckverlauf 
5 % 29 Be P—m=—0 U} &cos DE a Re I (6.28) 
“ und der Wandkontur Mas 
B Br: = Ye Agcoske Fo). nein: BT 
E. (Yw war wie alle y„-Abmessungen mit ö, dimensionslos gemacht.) Der Beitrag einer Wellenlänge 


‚ von Z bis L+ A der welligen Wand (Breite wieder B) wird 


27 
Wı(4; 2) = Böu|pdy 


Be EN TE er N ( 6.30) 
00, prZFE A si 
? = z Arsin p 
Wegen Aysing@= Abö, (siehe (6.20a)) und (6.25), sowie wegen ö,/L= Re/2(L) wird | 
e : 
RR N (6.31). 


2& 


lee DIET Beitrag der Welle von Zbis L-+ A zum Gesamtwiderstand ist W (A) W,(4; Z) +W) (};L) 
' und also nach (6.24) und (6.26), (6.31): 


WARST Wald.) ea 5 “ (6.32). 
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Die gegenüber der ebenen Platte sich ergebende Verminderung des Reibungs- 
widerstandes durch Oberflächenwelligkeit wird gerade durch den zusätzlich 
hinzutretenden Druckwiderstand kompensiert. 

Diese Aussage gilt in zwei Ordnungen bezüglich e. Wir können das Ergebnis sorgfältiger 


auchyso formulieren, daß der Widerstand W (A; L) der Welle sich von dem Widerstand Wr (A; Z) 
des entsprechenden ebenen Plattenstücks erst in 3. Ordnung bezüglich e unterscheiden kann: 


E- wie Wr UD= Ole) miin23. (6.33). 


$ 7. Zusammenfassung. 

Zur Untersuchung des Einflusses von Welligkeiten umströmter Wände auf den Verlauf der 
laminaren Grenzschicht wird eine Theorie entwickelt, die sich im Sinne einer Störungstheorie der 
1. Ordnung auf schwache (periodische) Welligkeiten bezieht und die ferner von Vereinfachungen 
Gebrauch macht, die sich mathematisch für kleine Wellenlängen A (klein relativ zur überströmten 


ET I NEE EEE ee rn 
Mer e ar Kr AT RER NER SPORE er ne 
, x Ray. Fe Te 
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= Wandbogenlänge L) erge 
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j Unter den genannten Voraussetzungen und soweit die Grenzschichttheorie überhaupt an- 
 wendbar bleibt, ergibt sich das folgende Verhalten. Die Welligkeit der Wand wirkt sich sehr stark 
auf den Verlauf der Grenzschichtströmung aus. Die Störungen gegenüber dem Fall der ebenen 


Platte sind in nächster Wandnähe relativ am stärksten. Am äußeren Rande der Grenzschicht ist 
- jedoch die Stromlinienwelligkeit gegenüber der Wandwelligkeit bereits sehr stark abgemildert. 


ie Be, 
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Diese Merkmale traten bereits in den von Quickund Schröder[9]berechneten numerischn 


Beispielen zutage. Die Dicke ö, der wandnahen „Störungsgrenzschicht‘, in welcher durch die 
Störung verursachte Reibungskräfte wirksam werden, ist klein gegenüber der Dicke ö der unge- . 


störten Grenzschicht, und zwar gilt das einfache Gesetz, daß - j 


d, 20 JAL 


ist. Wenn nur genügend Wellenlängen überströmt werden, kommt es zur Ablösung der Grenz- - 


schicht (wie auch schon Quiek und Schröder an dem Verhalten der von ihnen berechneten 

numerischen Beispiele demonstrierten). Zunächst erfolgt lediglich in allernächster Wandnähe 
innerhalb der kleinen Störungsgrenzschicht Ablösung mit Wiederanlegen, wobei im abgeschlosse- 

nen Wirbel des Rückströmungsgebietes praktisch Totwasser herrscht. Für das Einsetzen der Ab- 

lösung lassen sich einfache Gesetze formulieren. Sie interessieren im Hinblick auf-den laminar- 

turbulenten Umschlag. Das angegebene Ablösungskriterium beansprucht als Folgerung aus einer 

Störungstheorie der 1. Ordnung nur größenordnungsmäßige Richtigkeit. Es wird durch die 

Quick-Schröderschen Beispiele befriedigend bestätigt. 


Eine erste Ausdehnung der Untersuchung auf eine Theorie der 2. Ordnung im Störungs- 
parameter e (Amplitude der relativen äußeren Geschwindigkeitsschwankung) zeigt, daß der Rei- 
“ bungswiderstand der welligen Platte kleiner wird als derjenige der ebenen Platte, in Übereinstim- 
mung mit den Feststellungen an numerischen Beispielen von Quiekund Schröder. Estritt 
jedoch ein Druckwiderstand hinzu, der diese Widerstandsverminderung im Rahmen der Theorie 
2. Ordnung gerade wieder kompensiert. Eine Änderung des Gesamtwiderstandes durch Welligkeit 
der Oberfläche kann frühestens in der 3. Ordnung in & entstehen. 
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Die Stabilitätsfrage bei rückgekoppelten Verstärkern. I) 
| Von N . Joachim Lehmann in Dresden. - Re rt | E Br f 


3 Im vorliegenden ersten Teil der Arbeit werden die Grundlagen zur rechnerischen und zeichnerischen 
‚ (experimentellen) Stabilitätsuntersuchung linearer, rückgekoppelter Verstärker von einem einheitlichen - 
SEE mathematisch-physikalischen Standpunkt betrachtet und in eine für die Anwendung geeignete Form gebracht; 
der später folgende zweite Teil bringt dazu Beispiele (Vierpoltheorie) und Ergänzungen. 


‚Inthe present first part of the report, the foundations ofthe numerical and graphic, (experimental) methods 9 

=... of inwestigating the stability of the linear amplifyer the output of which is connected to the input are consi- 2 

en dered from a unitary mathematical-physical point of view, and transformed in such a way that they may ! 
be easily applied. A second part will follow containing ewamples (theory of quadrupoles) and supplements. 


Dans la premiere partie du prösent essai om etablit les bases de l’etude numerique et graphigue (experi- 

mentale) de la stabilite des amplificateurs lineaires & reaction, consideres du point de vue mathematique et 

Br physique, en meltant la theorie sus une forme appropriee aux applications pratiques. La deuaitme partie Sk 

de l’essai, qui doit pasaitre plus tard, en donnera quelgues ewemples (theorie des quadripöles), ainsi que les Be 
complements. a 


B mepsof uacru nammof padorTsI PaccMaATpuBamTcH C (PH3NKO-MATEMATnuecKofi TOYEH a 

- 3peHHA TEOPETHYeCKUe OCHOBAHHA UYHCHEHHOTO H TPahnyeckoro (3ECHEPHMEHTANBHOTO) HCCHE- Kr 

AOBAHNA. YCTOÜYUBOCTNU JIMHeÜHBIX "YcusuTeneli .C OÖPATHOH CBASbBE, Ipmuyem HTH Teope- 

THYeCKHe OCHOBAHUA IIPHBONATCH B YNOÖHYIO MIA npmMmeueHnun dopmy. Bo sropoü yacru Be. 
< pa6oTst ÖynyT naHsI ITPHMepsBr MH AONOIIHEHHA. 


= Einleitung. 

4 2.0 Die Frage, ob ein Verstärkergerät stabil weiterarbeiten kann, wenn es rückgekoppelt wird, 
-- Mist von verschiedenen Seiten theoretisch behandelt worden. Die besonders für die elektrische 
@ Nachrichtentechnik grundlegenden Arbeiten stammen von Strecker?) und Nyquist°), 


wobei allerdings gegen die Nyquistsche Auffassung des Selbsterregungsvorganges schwer- 
s wiegende Einwände vorgebracht wurden ®). Strecker behandelt die Aufgabe sehr allgemein, 
2 trotzdem fehlt aber eine einfache anschauliche Erklärung verschiedener scheinbarer Paradoxa, 
4 wie er gelegentlich selbst betont ®); außerdem wird infolge der vorwiegend anschaulichen Betrach- 
tungsweise übersehen, daß die benutzten Vierpolkenngrößen nur unter Nebenbedingungen zur 
Stabilitätsuntersuchung ausreichen 5). In der folgenden Arbeit soll die gesamte Theorie der Selbst- 
“erregung nochmals in einfachster Fassung dargestellt werden ®); durch Beschränkung der Unter- 
suchungen auf eine einzelne physikalische Kenngröße des Verstärkers wird das Verständnis der 
Vorgänge erleichtert und eine gewisse Vollständigkeit bei der Betrachtung der in der Praxis vor- 
kommenden Möglichkeiten erreicht. Es zeigt sich hierbei sehr deutlich, daß die Rückkopplung 
keineswegs mehr als Verstärkungs- sondern als Eigenschwingungsproblem betrachtet werden muß, 
und daß die Nichtbeachtung dieser Tatsache der Grund für alle scheinbaren Paradoxa ist. Bei 
der Betrachtung des Vierpols ?) werden mehrere für die Stabilitätsuntersuchung geeigneten Kenn- 
größen angegeben. ? 


F Teil L: 
Die eindimensionale Theorie der Rückkopplung.’) 


1. Der Verstärker. 


Es wird hier ein linear arbeitender Verstärker (Übertrager) Verstärker _ Belaslung5syslem 
‚„„V“ vorausgesetzt, der am Ausgang durch ein angekoppeltes System 
„B“ belastet sein kann (kurz: Belastungszustand B des Verstär- BAM einge age ee 
kers; in Bild 1 ist die zugehörige Schaltung als Meßschaltung be- are ee Dee 
zeichnet). Das System „B“ soll dabei in gleicher Art wie der Ver- Verbindung (172; 172) = Rück- 


= x 5 kopplungsschaltung. 
stärker ‚„‚V‘“ linear arbeiten. 


1) Mitteilung des Lehrstuhls f. angew. Mathematik (Prof. Dr. Fr. A. Willers) der T.H. Dresden. 

2) Anschauliche Theorie der Selbsterregug von Felix Strecker, Hirzel, Leipzig 1945. Darin ist ein un- 
veröffentlichter Aufsatz aus dem Jahre 1931 (am 28.4.1931 bei der Zeitschr. ‚Elektr. Nachr. Techn.“ eingegangen) abge- 

- druckt, der den Prioritätsanspruch des Verfassers begründet. (Im folgenden mit Strecker zitiert.) (Das Buch ist 

noch nicht erschienen.) 

») H. Nyquist, Regeneration Theory, Technikal Publication. Bell Telephone System 1932. 

“) Strecker, Abschnitt O3. _ 

5) Man vergleiche hierzu TeilII, Abschnitt 2 u. 3 dieser Arbeit. 

6) Die vorliegende Arbeit war bereits im wesentlichen fertiggestellt, bevor Einsicht in einige Fahnenkorrekturbogen 
(bis 8.64) des Streekerschen Buches genommen werden konnte. ; t 

?) Obwohl die Theorie im folgenden ganz allgemein physikalisch-mathematisch ohne Bezugnahme auf spezielle An- 
wendungsgebiete entwickelt werden soll, werden viele Bezeichnungen aus der elektrischen Nachrichtentechnik benutzt: 
z.B. wird eine Zeitfunktion u(f) in Anlehnung an den Röhrenspannungsverstärker mit „‚Spannung‘‘ bezeichnet. 


ve 
u 
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24 Lehmann, Die Stabilitätsfrage bei rückgekoppelten Verstärkern 4 
” er BR). 
Der Zusammenhang zwischen der angelegten Eingangsspannung ı, (?), die als Funktionder 
Zeit t alle später benötigten Differentiierbarkeitseigenschaften haben soll, und der Ausgangsspan- L 
nung %, (%,) = f (t) sei durch zwei lineare Differentialoperatoren | = 
. n d' m rg ; 
ne No mtc,+0.und M= DS D- mi D. 0 
p2 re 2, dt’ $; 
gegeben: 
Nu,=-Mm Dune sa & Dee Re 2 re 


Diese Differentialgleichung gilt für einen bestimmten Belastungszustand B; die ©,, D, und damit 
auch die Bauelemente des Verstärkers sollen zeitlich konstante Größen sein. Durch GI. (3) ist 
die Arbeitsweise des Verstärkers nur in bezug auf die Bewegungsgröße u (t) mathematisch gegeben ; 
das wird aber bereits genügen, um die charakteristischen Merkmale der Rückkopplung und das 
Wesen der scheinbaren Paradoxa zu erkennen. Beispiele und die Anwendung bzw. Erweiterung 
der Ergebnisse auf ‚‚mehrdimensionale Verstärker‘ (Vierpole), bei denen gleichzeitig mehrere 
Zeitfunktionen « (f), v(t), ... betrachtet werden müssen, werden im Teil II durchgeführt. 
Sofern die Betriebseigenschaften des mit Gl.(3) gegebenen Übertragers durch Messungen fest- 
gestellt werden sollen, wird er im allgemeinen ‚‚stabil“ arbeiten müssen, d.h. «, muß nach Ab- 
schalten von %, wieder auf Null zurückgehen: 


lim w, (u, () = 0) Le (4) 
t—> » 


i 
‘ 
| 


Ist aber diese Bedingung erfüllt, so haben die zum Operator M gehörigen charakteristischen 
Zahlen nur negative Realteile. Die Lösung von (3) lautet für w, (ft) = 0 für t>0: 


a 
u, (u, = 0) = Be Pe ER TE ds (5), 
e=1 
dabei sind die charakteristischen Zahlen 2„,.= Öm,oe Fiom,.(e=1,..,a <m) die «a vonein- 
ander verschiedenen Wurzeln der charakteristischen Gleichung: 


% DES LD, = TED, 2.— a = MW ld) (6), 


v=0 e=1 


fi gibt ihre Vielfachheit an (m= ) 1). Physikalisch stellen die z„, , das im allgemeinen 


e=1 
komplexe Anklingmaß der Eigenschwingungen an der Ausgangsseite des Verstärkers dar. 
Die in GI. (5) auftretenden Ausdrücke P;(—-V sind Polynome (#,—1)-ten Grades in t: 
un—1 


4 
PR — > ee (7). 


o=0 


Die m Konstanten p,,. werden durch den Anfangszustand von %, zur Zeit = 0 vollständig 
bestimmt. 

Damit die Stabilitätsbedingung GI. (4) für jeden beliebigen Anfangszustand erfüllt.ist, folgt 
aus Gl. (5) offensichtlich die Behauptung 


N, 01 na RE (8). 


Sobald der betrachtete Verstärker ‚stabil‘ arbeitet, kann man leicht seine Übersetzung Ü (io) 
messen. Diese Kenngröße wird besonders in der Elektrotechnik benutzt und ist als Verhältnis 
der Ausgangsspannung im eingeschwungenen Zustand (d.h., «,=f(t) wird nur bei t— oobe- 
trachtet) zu einer sinusförmigen Eingangsspannung mit der Kreisfrequenz o definiert. D.h., wenn 
man in komplexer Schreibweise u, = E et®t setzt: 


U Part (Bei® r) 


U (io) == naar ar a EEE EEE (9). 


Das hier auftretende Partikularintegral %,part (Beiet) der GI. (3) wird unten bei Gl. (11) 
genauer betrachtet. 

Die durch Gl. (9) gegebene Kennfunktion für den Übertrager genügt meist vollständig’ zur 
Charakterisierung seiner Verstärkereigenschaften, sogar die Stabilitätsfrage des rückgekoppelten 


Den °) Vgl. z.B. L, Bieberbach, Theorie der Differentialgleichungen. 3. Aufl., S. 332/336. Berlin: Julius Springer, 


gelegt . N zwar technisch Schwierigkeiten, mathematisch setzt, 1 man einfach Ar 


. h 


Bau Fa VE isn, U ae ee we Kai ee 


3 23 yet 
2m, o 


den BES WÜRSP ES En zu Gl. (9 bestimmen: 


Uq Part (E e? e*‘) 


Tas Wert ee ee 


eu Malsancch stellt das einfach die analytische Fortsetzung von Ü (io) in die z-Ebene hinein: 


E:- ” . dar. ‚Für. Ya Part findet man aus Gl. (3) mit «,—= Be*! und dem Ansatz „= Be*: 
R- DE (Bei) "ID Ka Be a a a He ee, nn 49), 
E: wenn entsprechend GI. (6) 
= - n b 7 
2 ra 1 BO ae OB ee (14) 
4 Fe | . E% =0 e=1 e 
gesetzt wird. Damit wird die nach Gl. (12) definierte Übersetzung Ü (2): 
E ER 
U ; ; 
On (15) 
2 Diese gebrochen rationale Funktion in z hat höchstens in den Nullstellen von IM (z) Pole; 
En. physikalisch heißt das, nur für Schwingungen mit einem Anklingmaß 2 — 2m, . einer Eigenschwin- 


gung der Übertragerausgangsseite kann die Übersetzung über alle Grenzen wachsen. Im Gegensatz 
zu der mit Gl. (12)—(15) eingeführten (verallgemeinerten) Übersetzung Ü (z) bezeichnen wir die 
nach Gl. (9) gemessene Kenngröße Ü (io) als w-Übersetzung. 


‘Nach diesen Betrachtungen am Verstärker führen wir einen Rückkopplungsweg ein. 


NT Tee) Khan 


4 

; RE | B | - 2. Die Rückkopplung. 

Es wird jetzt der Eingang des unbelasteten Verstärkers mit dem Ausgang verbunden, so daß 

| zwangsweise %, = U, = U wird (hierzu Bild 1) °). Zur mathematischen Formulierung dieser Rück- 

kopplung gehen wir von der Differentialgleichung (3) aus. Dabei müssen wir aber den Übertrager 

beieiner Belastung R betrachten, die der zwangsläufigen Belastung beim Zusammenschalten von 

7 -  Ein- und Ausgang bei der Rückkopplung entspricht. Wenn dieser Belastungszustand durch einen 

Index R angedeutet wird, beschreibt die homogene lineare Differentialgleichung 


NUN EN BT I NE MET re in 


das Verhalten des rückgekoppelten Verstärkers1®) — der rückgekoppelte Apparat führt nur ‚‚Be- 
wegungen“ u(t) aus, die sich aus gewissen Eigenschwingungen zusammensetzen lassen: 


= I’ PWe=D ern.«. a a FON 


e=1 


®) Überlagert man am Eingang noch eine Zwangsbewegung uz(t), so daß u, = uy+ u, wird, bekommt man an Stelle 
der homogenen Differentialgleichung (16) eine inhomogene: N RUE M RU —N Rlde die zugehörige Übersetzung ist 


entsprechend Gl. (12) u. (15): a - NR (MEI —NR@): 
16) Hierzu liefert die ‚‚Vierpoltheorie‘‘ wichtige Beispiele, man vgl. II, 2 u. 3. 


ates kan t entschieden en Vesnduie der Zusaneeane beim rück- 
pelten Verstärker kommt man aber nur durch eine Erweiterung des Übersetzungsbegriffes. je x 
ierzu denken uns an den Verstärkereingang eine sinusförmige an- oder abklingende Spannung Ro: 


Er = = 9 5 io LE komplexem Anklingmaß. Am Von „mißt‘“‘ man nur den En i 
er nungsanteil mit dem Anklingmaß der Eingangsspannung: %, part (Bet) = (2); Eigenschwingun- 
gen des Apparates an der Ausgangsseite, die nach Gl. (6) das Anklingmaß z,,, haben, bleiben 
damit unberücksichtigt. Rau 2 - zufällig mit einem ms R überein, so wird UyPart en Ad %) dur an 


Ysrart (Beim,ei) — im UyPart (Beirat EEE a) Kr; = 


3 efinlert, Damit kann man eine verallgemeinerte Übe rsetzung Ül (2) bei an- und abkdingen- 
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 ristischen Gleichung des rückgekoppelten Verstärkess 
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-ß,ist wieder ihre Vielfachheit. Die in den entsprechend Gl. (7) definierten 'olyno: 
«Pe» der Gl. (17) auftretenden Konstanten P,,« ‚werden durch den Anfangszustand. 


lediglich durch seine r verschiedenen Eigenschwingungen unter Berücksichtigung von deren i 
‘den ohne Änderung der Schaltung von außen nicht mehr eingegriffen werden 


dert aber nichts daran. Die Übersetzung Üz (2) — d.h. die Übersetzung eines offenen Verstärkers 
bei der Belastung R — hat nur noch mittelbare Bedeutung: Üz (z) = 1 entscheidet über die Sta- 
- bilitätseigenschaft des rückgekoppelten Gerätes, während ü= Üz/l—Ü, dessen Übersetzung 


+ 
+ 
Er Nr 


.,=0 v=0 fiz 


= 


Po 


en 0 .festgelegt. 3 A , N Rt, 
An dieser Stelle erkennt man deutlich, wie das Verhalten des rückgekoppelten Verst 


fachheit bestimmt ist: derikt man sich z. B. zur Zeit = 0 am Eingang des rückgekoppelten V: 
stärkerseinen bestimmten Anfangszustand «, (0), « (0), ... vorgegeben, so werden dadurch die] Bon 
stanten ?,,. im Ansatz Gl. (17) bestimmt. Für t>0 läuft ein Eigenschwingungsvorgang ab, n 


A 


kann; es darf — auch nicht im Gedankenexperiment — der Übertrager-Ein- und Aus- BR 
gang nicht mehr getrennt beobachtet werden. Der Einfluß einer am Eingang zusätzlich auf- 
gezwungenen Bewegung 4z (t) überlagert sich diesem Vorgang einfach linear (s. Fußnote ®), än- 
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nu 
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ist (s. Fußnote ®). Nichtbeachtung dieses Tatbestandes ist der Grund für alle scheinbaren Para- 
Aoxa-y, x See 


Eh 
| 3. Die Entscheidung der Stabilitätsfrage. re Fa 
Die Schwingung % (t) nach Gl. (17) — und damit der rückgekoppelte Verstärker — wird 
instabil, sobald eine der Wurzeln zz,. der zugehörigen charakteristischen Gl. (18) einen positiven 
Realteil hat; bei nur negativen Realteilen klingen alle Teilschwingungen ab. Ein stationärer Vor- 
gang ist nur dann möglich, wenn wenigstens eine Zahl zz,, rein imaginär und einfach ist und alle 
anderen negative Realteile haben. Bei der Untersuchung der Stabilitätsfrage sollen zwei Fälle 


unterschieden werden: 5 & 
A.Es ist die Differentialgleichung (l6)bekannt. Be, 
In diesem Falle muß die ganze rationale Funktion Nx (2) — Mz (z) nach Gl. (18) aufgestellt r 


und auf Wurzeln mit positivem Realteilgeprüft werden. Dazu stehen die inhaltlich gleichwertigen 
Kriterien von Routh 2) und Hurwitz!?) und ein anderes von Leonhard!) zur Ver- . 
fügung. Bis Gleichungen 6. Grades hat dabei Bilharz ») eine graphische Methode angegeben, 
die eine einfache Handhabung des ersten Kriteriums ermöglicht. In ihrem Anwendungsgebiet 
erlaubt die graphische Methode von Bilharz die schnellste Stabilitätsuntersuchung, allgemein 
stellt wohl das leicht zu schematisierende Routhsche Verfahren *) die einfachste Prüfungs- 
methode dar 17). - ; 


‘ 


B. Es wurde die @-Übersetzung Üz(io) experimentell aufgenommen. 


Hierzu betrachten wir den Zusammenhang zwischen der Hauptgleichung (18) und der Über- 
setzung U, (z) genauer. Division der Gl. (18) durch M; (z) ergibt: 


GA TI- RN SFEEIRER 


Diese Gleichung liefert zumindest alle Wurzeln zz, , von Gl. (18), die nicht gleichzeitig Null- 
stellen von WM; (2) sind, in der richtigen Vielfachheit. Und das genügt zur Untersuchung von ° 
Stabilitätsaufgaben aus der Praxis wohl immer, besonders, da hier nur Lösungen mit positivem 
Realteilinteressieren und nach Gl. (8) die Wurzeln von M (z) bei einem „‚stabilen‘‘ offenen Verstär- 
ker immer negative Realteile haben. Da die Wurzeln der ganzen rationalen Funktion Nz (2) — Mp (2) 
und damit auch die von Gl. (19) in der z-Ebene alle innerhalb eines Kreises um den Nullpunkt mit 
einem endlichen Radius $ liegen 18), so ist die Suche nach Wurzeln mit positivem Realteil gleich- 


11) Man vgl. hierzu auch das Beispiel in II, 1. 

=) E. J. Routh, A treatise on the stability of a given state of Motion. London 1877, 

») A. Hurwitz, Math. Ann. Bd.46 (1895), 8.307—311. Das Hurwitzsche Kriterium wurde von H. Bil- 
harz, Z. angew. Math. Mech. Bd.24 (1944), S.77—82, für Gleichungen mit komplexen Koeffizienten erweitert. 

“) A. Leonhard, Arch. Elektrotechn. Bd. 38 (1944), S. 17—28. 

‘) H. Bilharz, Z. angew. Math. Mech. Bd. 21 (1941), S. 96—102. . une 

*) Fr. A. Willers, Prakt. Analysis. 2. Aufl., $30. Berlin: de Gruyter. Rechenformulare hierzu liegen z.B. 
im Math. Seminar der T.H. Dresden. 

ı?) Ein Beispiel hierzu wird in II, 4 durchgerechnet. 

"®) Man vgl. hierzu und im folgenden z.B. K. Knopp, Funktionentheorie. Sammlung Göschen, Bd. 668. 6. Aufl. 
1398. Berlin: de Gruyter 1944. 
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sliegen, welches 
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euten 


d mit der Frage, ob in der z-Ebene Lösungen z;, , innerhalb eines Gebiete 


... durch eine Kurve J entsprechend Bild 2 begrenzt wird. (Der Weg J setzt sich allgemein aus dem 


- Wurzeln 25, = io@r,: oder Pole 2, 
umrandeten Gebiete ausgeschlossen werden 19), 


Stück von —S bis +8 der imaginären Achse und einem Halbkreis um den Nullpunkt in der 


= 
Ma 


ö> 0-Halbebene mit dem Radius 8 zusammen. Sollten dabei zufällig auf der imaginären Achse 


pn, e liegen, so sollen diese durch einen Haken aus dem von J 


Bei der in der Praxis vorliegenden Aufgabe ist Ür (2) nur auf einem Teil der imaginären 


Achse der z-Ebene bekannt (Ü7 (io)); vorläufig sehen wir aber einmal die vollständige durch Üz (z) 
. vermittelte Abbildung der J-Kurve in eine Ü=X-+iY-Ebene als gegeben an. Da Uz() auf J 


überall regulär analytisch ist, wird J in der Ü-Ebene wieder als geschlossener Kurvenzug J’ mit 
einem bestimmten Richtungssinn abgebildet. Die Funktionstheorie liefert nun eine einfache Vor- 


schrift zur Prüfung, ob der Kurve J’ im Innern der J-Kurve Wurzeln 25, „entsprechen oder nicht?°). 
E.. Es handelt sich dabei um eine bekannte Anwendung des Residuensatzes: Es sei f (z) in @ ein- 


_ 


Bild 2. Die z-Ebene mit der Bild 3. Die Ü-Ebene mit der Abbildung der 
J-Kurve. J-Kurve und einem Teil der ö-Konst. bzw. w- 
Konst.-Linien. (Es ist Ü (2) = (1—i1,5 + 2)! 
’ gewählt. e 


deutig (evtl. auf einer mehrblättrigen Riemannschen Fläche) und regulär, C’ist ein in G liegen-- 
der geschlossener doppelpunktfreier Weg mit dem J entsprechenden Richtungssinn. Ist dann 


- f(z) auf C von Null verschieden und liegen im Innengebiet von © höchstens endlich viele singu- 


läre Stellen, die sämtliche Pole sind, so ist - 


(©) : 
: a a (20), 
£ 2,8). fi2) 
wobei N die Zahl der Nullstellen, P die Zahl der Pole von f (z) innerhalb O angibt; jede Stelle wird 
dabei ihrer Vielfachheit’entsprechend gezählt. Setzt man nun f(z) = R (x, y) ei?» und damit 
f'(z) dae=eiPdR +iReiPdD, wobei R (x, y) in einer f-Ebene der Betrag, ® (x, y) das Argu- 
ment von f (2) ist, so wird aus Gl. (20) 
i Far 1 er 1 ie 
ER == 5 
j i Ir ler fee je nen 
Das erste Integral fällt hierin weg, da nach Voraussetzung immer R (x, y) reell und >0O 
bleibt und über einen geschlossenen Weg integriert wird. Der Satz Gl. (21) wird jetzt auf die 
Funktion f(2)= Üz (z)— 1 und den Integrationsweg J bzw. sein Bild J’ in der Ü-Ebene ange- 
wendet. Da f (z) hierbei eine gebrochen rationale Funktion ist und J Nullstellen und Pole von 
j(z) vermeidet, sind alle Voraussetzungen erfüllt, und man erhält: 
(N) . 


1 
a ED Re ee ee 
P-n=—| (22) 


In der Ü-Ebene ist hierbei ® wegen f (z) = Re’? = Üz(z)— 1 der Winkel eines Vektors 
im Punkte Eins. (Bild 3.) 


12) Im Hinblick auf die praktische Anwendung wurde der hiermit für J gegebene Durchlaufungssinn als positiv an- 


gesehen." 
2) K. Knopp, a.a. 0O.,.S.127—131. 
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_ schließt nl bekommt dann das Er von Nyquis 


keit, daß Pole auf der imaginären Achse liegen, etwas genauer betrachtet). 


(Beim stabilen rückgekoppelten Verstärker muß N = 0 sein!) 
Man zeichne in der Ü-Ebene vom Punkt Ü= 1 aus einen Vektor R,e 
ne Punkte der J’-Ortskurve. Läßt man dann den Endpunkt des Vektors ein 
J’-Ortskurve herumwandern, so ändert sich der Betrag des Winkels ® gegenüber m 
winkel ®, dann und nur dann um Vielfache von 2 x, wenn das System instabil ist. 
Der Fall, daß nicht P= 0 vorausgesetzt werden kann, kommt bei „mehrdimen: 
Verstärkern“ vor und wird daher beim Beispiel des Dreipols, i in II, 3, ‚zusammen or Ei 1 


' Der praktischen Anwendung des mit Gl. (22) bewiesenen Unitänfieriterkunie 
Falle, daß P = 0 und das Fehlen von Polen auf der ar Achse vorausgesetzt ‚kan 
noch viele Schwierigkeiten entgegen, denn es wird ja nur Üz(io) in einem begrenzten o-Interval 1% 
gemessen. In der S. Ebene ist also nur ein Teil der Abbildung der imaginären Achse der z-Ebene 
bekannt; die J’-Kurve ist nicht vollständig gegeben. Um hier weiterzukommen, beachten wir 
das asymptotische Verhalten und die Konformität der durch Üz(z) gelieferten Abbildung; dabei 
wird auf der ÜR(iw)-Kurve immer die &-Skala als bekannt angenommen. Zudem soll gleich 
berücksichtigt werden, daß bei den allermeisten in der Praxis auftretenden Aufgaben die in Pr 
auftretenden Konstanten reell sind; in der Elektrotechnik ‘handelt es sich dabei um 

wie Widerstände, Kapazitäten usw. In diesem Fall vereinfacht sich die Untersuchung be- 
trächtlich, obwohl auch der allgemeinere Fall ganz entsprechend behandelt werden kann. Bei 
reellen Konstanten ist die Übersetzung eine reelle Funktion der komplexen Variablen z, sie liefert 
also zu konjugiert komplexen Punkten z, z konjugiert komplexe Werte Ü, Ü. Man könnte sich De 
daher bei der Stabilitätsuntersuchung darauf beschränken festzustellen, ob im ersten Fre, x 
der z-Ebene Lösungen der Hauptgleichung liegen. Es ist jedoch genau so einfach, die theoretischen _ 
Überlegungen mit der J-Kurve durchzuführen und erst später zu berücksichtigen, daß man auch 
mit deren imersten Quadranten liegenden Teil auskommt; man kann den WegJ zerlegen: JIaIhJ, - 
wobei J, das Geradenstück von Oi und den Viertelkreis von 38 vertritt. JJistdrandr 
EV re gespiegelte, also zu J, konjugiert komplexe Weg; das negative Zeichen davor zeigt an, 
daß J, mit negativem Richtungssinn durchlaufen werden muß, da J, bei der Spiegelung von Jı 


_ den Richtungssinn von 0>— iS usw. zugeordnet bekommt. Da Üz (z) hier eine reelle Abbildungs- Be 


funktion der komplexen Variablen z ist, wird auch J’= J/—J' und damit wegen der allgemein 2 
(4) nn 


‚gültigen Beziehung [d® = — F. d®: 


Ke  ET  a E (07) 67) \ 


fa9=[a9= [ao [ao=2[a8.......... 23). 


Diese Gleichung ermöglicht die endgültige Formulierung der Stabilitätskriterien ohne Be- 
zugnahme auf Üz (iw) bei negativen Frequenzen ». Wenn die w-Übersetzung als gegeben ange- 
sehen wird, fehlt an der J’-Kurve noch das Bild des Halbkreises von J, dessen Radius S auch un- 
bekannt ist. Es wird der Versuch unternommen, diese Lücke zu schließen. Dabei unterscheidet - 
man zweckmäßig zwei Fälle: 


a)m, der Grad von Wx(), istgrößerodergleichn,dem Gradvon a 
In diesem Falle verhält sich Üz (2) im Unendlichkeitspunkt regulär, es ist lim Üz (2) = 


wobei für m>n Ü„=0, sonst ein im Endlichen liegender Punkt ist. Das Bild des Halb- 
kreises der J-Kurve schrumpft dabei mit wachsendem Radius Sin den Punkt Ü, hinein, so daß 
die J-Kurve in der Ü-Ebene mit $— o bei Kenntnis von Ü, (io) geschlossen werden kann. Be- 
RIchdet man das in der Ü-Ebene über den Weg Ü, (io) von = 0 x geführte Integral ; d® 


it [d®, so wird bei reeller Abbildungsfunktion Üz (2) mit S> oo (im Weg J) aus Gl. (22) 
ar: (25): 


(0) 
PEN EN 
IT 


woraus z.B. im Falle des ‚‚stabilen eindimensionalen“ Verstärkers wegen P = 0 sofort die Sta- 
bilitätseigenschaft des Systems entnommen werden kann. Ob der Fall m>n vorliegt, erkennt 


»!) Obwohl ein offener eindimensionaler Übertrager mit Polen auf der imaginären Achse in der Pr. 
axis meßt 
kaum stabil arbeitet, vgl. man hierzu z.B. das Beispiel nach Bild 10 u. 11 im Artikel von A. Leonhard,a. . 


ÖL) Urn) Bl, dh, Grm... 
wo - asymptotische Gleichheit für lz1- — o angibt. Für genügend große. |z| gilt also eine Näherung = 
‚in der Form: 

k ie REN ET Re EEE (26). 
Für das Bild des rn der J-Kurve für sehr großen Radius 8 folgt mt 2— Sei?, 
"wobei p die Werte Von. — > us 5 durchläuft; 
Ün (Sein) IE N DA N NE N re Er Let (27), 


| Q 


F 
f-- 


ı @ bei schnell enger werdender -Skala dem Punkt Ü „ nähern, der bei m>n im Nullpunkt, 


bei m = n und reeller "Funktion ÜR (2) irgendwo auf der reellen Achse liegt. Restlose Klarheit 
bringt dann die analytische Fortsetzung der Ü, (i@)-Kurve zyerelatn auf > ©) 2) 23), Meist ee 
gt ein roher graphischer Überblick. aut 


b) Es ist meklheineralsn * ir he 


Ür (z) hat damit im Unendlichkeitspunkt der rrbene einen Pol (n-m)- -ter” Dr der 


durch den Halbkreis der J-Kurve aus dem Integrationsgebiet ausgeschlossen werden muß. 


„Hier zerlegt man Üz(z) in eine rationale Funktion Are) vom Grade (n— m) = 
und eine echt gebrochene Funktion ® (z): 


x und 7 die Werte A) rc bis a, ; annimmt. Da A sehr groß ist, wird auch #> ®, 


2 2 
und man bekommt für S—oo schließlich mit Gl. (23) und Bezeichnungen nach (Gl. 24) aus 
Gl. (22): 


Age m)Z 


(n n—m)z 


| 2 j i 
rn... »[a0+ | aw en E a 


Um diese Beziehung zur Stabilitätsuntersuchung verwenden zu können, braucht man nur 
die Größe (n— m), die sich aber meist leicht aus dem Verlauf von ÜR (io) im Bereiche großer ® 
entnehmen läßt. Aus Gl. (25) wird mit 2=io: 


ÜR(io) > Kim nm) 4 K,in—m—1) an—m—1) +..= 4etr, 
5 d.h. |y|=0 oder r bei geradem, 5 bei ungeradem (n— m). Das ergibt die Regel: 
Ay Ar gerade : Ä f reellen 
Es ist (n— m) es ” wenn die w-Übersetzung mit wachsendem & der imaginären 


Achse zustrebt. Die genaue Größe von (n—m) ist dann einfach aus der w-Skala (die genähert 
am) darstellt) abzuschätzen, zumal (a—m) in der Praxis meist klein ist. Auch hier kann 
gelegentlich die analytische Fortsetzung mit Vorteil benutzt werden. 


Wurde bisher immer das Umlaufkriterium Gl. (22) zu Grunde gelegt, das zu einfachen Sta- 
bilitätskriterien führt (Gl. (24) und (28)), so treten doch Schwierigkeiten auf, wenn die Zahl der 
- Pole innerhalb J unbekannt ist. (Man vgl. II, 2). Es bleibt dann auf jeden Fall die analytische 
_ Fortsetzung der Ü,(iw)-Kurve in die z-Ebene hinein zur Untersuchung geeignet. Man entwirft 
dabeiin der Ü-Ebene (unter Berücksichtigung der Konformität der Abbildung) von Kurventeilen 
von Üz(iow) ausgehend, denen der kritische Punkt Ü= 1 rechts benachbart ist, das Bild eines 
Maschensystems von ö = konst., ® = konst.-Linien der ’z-Ebene (vgl. Bild2 u.3). Damit wird 
nicht nur festgestellt, ob dem Punkt Ü=1 in der z-Ebene Werte zz, , mit positivem Realteil 
entsprechen, sondern 25, , kann zahlenmäßig angegeben werden ®). 


22) Fr. A. Willers, Graphische Integration. Sammlung Göschen, Bd. 801. S.104—131. Berlin und Leipzig: 
- de Gruyter 1920. Weitere Literatur in Encyklopädie der Mathematischen Wissenschaften Bd.II, T1.III, 1. Hälfte im 
Artikel von ©. Runge u. Fr. A. Willers über numerische und graphische Quadratur und Integration gewöhnlicher 
und partieller Differentialgleichungen. 

2») Man vgl. hierzu z.B. Bild 3. 

2) Das Leonhardsche Stabilitätskriterium stellt eine Spezialisierung dieser allgemeinen Überlegungen auf 


ganze rationale Funktionen dar. 


25) Vgl. Fußnote 22 auf S.6. 


n n en len bei oberflächlic er Belang A, Üz (io)-Kurve; sie muß sich mit wachsen- 


(Bil 


Über die Machsche V-Ausbreitung. 


» r J 


_ Die V-Ausbreitung wurde von E.Mach!) bei der 


Durchdringung zweier Kugelknallwellen beobachtet, 
die gleichzeitig an zwei verschiedenen Orten A, B 
d 1) erzeugt sind. Zunächst ist die Durchc 


_ Bild 1. Durchdringung zweier Kugelwellen. 


so, wie es von Schallwellen her geläufig ist, wird aber 
ein gewisser Winkel 29 überschritten, so ändert sich 
das Bild, es entsteht nach den kinematographischen 
Aufnahmen von C.Cranz und H. Schardin?) im 
Punkt C eine neue Welle (Bild 2), die sich in den von 


Bild 2. Machsche ‚‚V‘‘-Ausbreitung. 


den Knallwellen noch nicht erfaßten Kielraum aus- 
breitet. Der Name dieser Erscheinung kommt von den 
Versuchen E.Machs mit berußten Glasplatten, auf 
denen zwei Knallquecksilberhäufchen zur Entzündung 
gebracht wurden. Im Gebiet der neu entstehenden 
We le ist der Ruß weggewischt, wodurch ein V-förmiges 
entrußtes Gebiet entsteht. 

Von C. Cranz und H. Schardin?) wurde das Kri- 
terium aufgestellt, die neue Welle entstehe, wenn der 
Druck im Überlagerungsbereich der Wellen sich 
schneller ausbreitet als mit der Geschwindigkeit des 
Schnittpunktes C (Bild 1). Hier soll dieses Kriterium 
gasdynamisch untersucht und für das Entstehen der 
V-Ausbreitung eine Erklärung gegeben werden. 

Wir betrachten zunächst die schiefe Reflexion eines 
Verdichtungsstoßes an einer ebenen Wand und machen 
den Vorgang zu einem stationären durch Überlagern 
der negativen Geschwindigkeit des Reflexionspunktes, 
Wir haben jetzt zwei schiefe Verdichtungsstöße. Über 
sie liefert das Stoßpolarendiagramm im Geschwindig- 
keitsplan alles Wissenswerte, da in unserem Problem 
Be andhelasmngen Aussagen über Geschwindigkeiten 
sind. 

Die Konstruktion des Geschwindigkeitsplans ist in 
Bild 3 durchgeführt. Es ist dabei von der Tatsache 


Bild 3. Schiefe Reflexion, stationäres Problem. 


Gebrauch gemacht, daß die kritische Schallgeschwin- 
digkeit «in allen Feldern, die von der Wand und den 
Verdichtungsstößen begrenzt werden, die gleiche ist. 


ı) E.MachundI.Wosyka: Wien. Ber. 72,1875; W.Rusciky 
ebenda 73, 1876; E. Mach und I. Sommer ebenda 75, 1877; 
B.Mach, 0.Tumlirz und C. Kögler ebenda 77, 1878; E.Mach 
ebenda 77, 1878. 

®) C.CranzundH.Schardin: Z.f: Phys.56 (1929) S. 147. 
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Welle und dem Einfallwinkel o zı 
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Bild 4. Geschwindigkeitsplan zur schiefen Reflexion. 


iz 


in dem mw; die Stoßpolare tangiert. Das bringt die 
bekannte Tatsache in Erinnerung, daß die Geschwindig- 
keit in einem schiefen Verdichtungsstoß nicht über 
einen bestimmten Winkel Omax hinaus abge 
kann. Wäre also der Einfallswinkel größer als der 

zu Ömax gehörige, so könnte die Geschwindigkeit inder 
reflektierten Welle nicht der Wand parallel abgelenkt 

werden. Infolgedessen besteht ein maximaler Einfalls- 

winkel omax. Er ist abhängig von der Intensität der 

einfallenden Welle. - 


Bild 5. Schiefe Reflexion beimaximalem Einfallswinkel. 
o= 40°, 0’ = 50°, 9 = 16°. 


Dieser stationäre Vorgang kann noch keine Erklärung 
für den nichtstationären Vorgang der V-Ausbreitung ° 
geben. Nimmt man aber die Wand nicht eben sondern 
gewölbt an, sohat man, wie im beobachteten Fall,einen ° 
sich stetig vergrößernden Einfallswinkel. Da nur die 
Entstehung der Y-Ausbreitung interessiert, genügt es, 
eine Wand mit einem kleinen Knick (Bild 6) anzu- 
nehmen. Der Verdichtungsstoß falle zunächst mit 
maximalem Winkel ein. Läuft dann der Reflexions- 


punkt über den Knick, so vergrößert sich der Einfalls- : 


lenktwerden 


- 


Zu ihrer Erklärung formulieren wir die Randbedin- 
gung der schiefen Reflexion anders, nämlich so, wie 
es beider Anwendung des Huygensschen Prinzips in 
der Akustik üblich ist. Die einfallende Welle bringt 
eine Geschwindigkeit U mit, die mit der Randbedingung 

nn -: nicht verträglich ist. Es muß also längs der Wand 
ü Weine Druckstörung einherlaufen, die die Normalkom- 
A| > Zu ‚ ponente von U. zum Verschwinden bringt. Die Ge-. 
'schwindigkeit w, — .D/sin omaz, mit der sie sich bewegt, 

ist größer als die Schallgeschwindigkeit hinter dem 
_ Verdichtungsstoß. Man wird so auf das bekannte 
Problem einer mit Überschallgeschwindigkeit bewegten 
Druckstörung geführt. Bei nicht zu großem Einfalls- 


 störung ausgesandten Wellen bilden eine ebene Front, 
die den reflektierten Verdichtungsstoß darstellt, ähnlich 
wie in der Akustik. Bei genügend großem Einfalls- 
winkel oa — und dieser besteht, sobald der Verdich- 
 tungsstoß den Knick der Wand überschreitet — ändern 
' sich die Verhältnisse, es bildet sich vor der Druck- 
störung ein senkrechter Verdichtungsstoß aus, wie es 
von.der Überschallströmung gegen einen Keil genügend 
großen Öffnungswinkels (Bild 6) bekannt ist. Damit is 


Bild 6. Entstehen der neuen Welle an der Kante einer Wand. 


die Entstehung des neuen Verdichtungsstoßes bei der 
Machschen ' V-Ausbreitung erklärt. Der neue Ver- 
: diehtungsstoß bildet mit dem einfallenden und dem 
| e. durch Reflexion entstandenen eine Gabel, denn das 
| Gas, das dem durch Reflexion entstandenen Verdich- 
tungsstoß nachströmt, hat gegenüber der Druckstörung 
Unterschallgesehwindigkeit, wie die Rechnung zeigt; 
CE hier kann sich also eine von der Druckstörung her- 
2% rührende Wellenfront nicht ausbilden. a 
Die Fortpflanzungsgeschwindigkeit der neuentstan- 
denen Stoßwelle ist i. A. nicht gleich der Geschwindig- 
> keit des Gabelpunktes. Bei schwachen Stoßwellen ist 
2 sie kleiner, bei starken Stößen größer. Daß trotzdem 
? die Gabel erhalten bleibt, steht dieser Feststellung, die 
sich auf stationäre Verhältnisse bezieht, nicht entgegen, 
denn die Gabel kann sich ändern. Bei schwachen Stoß- 
wellen ist der gerade Stoß, der sich mit der Geschwindig- 
keit des Gabelpunktes bewegt, zu stark, er sendet also 
Druckwellen aus, die die von früher her noch bestehende 
"reflektierte Welle verstärken und die Gabelform ändern 
(Bild 7), bei starken Stoßwellen sendet der gerade 
Stoß, der nun hinter sich einen zu kleinen Druck er- 
zeugt, wenn er sich mit der Geschwindigkeit des Gabel- 
punktes bewegt, Verdünnungswellen aus, die ebenfalls 
die Gabel ändern. Bei kleinen Druckdifferenzen sind 
dies akustische Wellen. Der Unterschied zwischen den 
beiden Fällen in Bild 7 kommt daher, daß im ersten 
ID Fall hinter dem geraden Stoß Unterschall- im zweiten 


Prof. Dr. 6. DOETSCH, Tabellen zur La- 
j place-Transformation und Anlei- 
'F tung zum Gebrauch, (Die Grundlehren der 
| mathematischen Wissenschaft in Einzeldarstellungen. 
- Bd.LIV.) VIII + 1858. mit 9 Abb. Berlin und Göt- 
| tingen 1947. Springer-Verlag. Preis geh. 18 M. 

Der Verfasser, dessen grundlegendes Buch über die 
Laplace-Transformationen vor 10 Jahren erschienen 
ist, gibt mit dem vorliegenden Bande dem Praktiker 
ein ‚wertvolles Hilfsmittel für die Anwendung dieser 
Transformationen in die Hand. Diese und ähnliche 


ulässige Maß, und es muß V-Aus- . Überschallströmung besteht. Wie die Nachrechnang * Ras” 


“ „winkel hat man den normalen Fall, alle von der Druck- 


=. ne —— 


zeigt, entspricht der Grenzfall der kritischen Strö- 
mungsgeschwindigkeit gerade dem Fall, in dem die 


“ 


Bild 7. Schiefe Reflexion. _ ö | 
Nichtstationäres Problem. > nr 


Geschwindigkeit der geraden Stoßwelle und die des 
Gabelpunktes übereinstimmen. Man sieht, daß das Kri- 
teriumvon Cranz und Schärdin durch das des maxi- 
malen Einfallswinkels zu ersetzen ist. 
Schließlich werde noch die Gasgeschwindigkeit U, 
hinter dem durch Reflexion entstandenen Verdich- 
tungsstoß und U® hinter dem neuentstandenen abge- 


Bild 8. Druckwelle vor einer mit Überschall-Geschwindigkeit 

bewegten Druckstörung. 
lesen für den Fall, daß die Geschwindigkeit des geraden 
Stoßes mit der des Gabelpunktes übereinstimmt. Nach 
Bild 4 ergibt sich UZ>U,. Das Machsche V erklärt 
sich jetzt so, daß hinter der neu entstandenen Stoßwelle 
der Ruß wegen der größeren Geschwindigkeit weg- 
gewischt wird. 


a) b) 
Bild 9. 


a) Schwache Stoßwellen, Überlagerung einer Druckwelle. 
b) Starke Stoßwellen, Überlagerung einer Verdünnungswelle. 


Der rotationssymmetrische Fall, in dem die einfallen- 
den und reflektierten Verdichtungsstöße Kegelflächen 
darstellen, unterscheidet sich, was die Entstehung der 
V-Ausbreitung und die Stärke des neu entstandenen 
Stoßes betrifft, in nichts von dem ebenen Fall, denn 
in den Grundgleichungen des Verdichtungsstoßes tritt 
die Krümmung der Stoßfront nicht auf. 

Aachen. F. Schultz-Grunow. 


Fe BUCHBESPRECHUNGEN 


Methoden haben sich in den letzten Jahrzehnten als 
Hilfsmittel zur Behandlung gewöhnlicher und parti- 
eller Differentiälgleichungen insbesondere bei den Elek- 
trotechnikern immer mehr eingebürgert, nicht zum 
wenigsten, weil sie Probleme, die auf partielle Diffe- 
rentialgleichungen führen, übersichtlich und nach ein- 
heitlicher Methode zu behandeln gestatten, während 
sie nach anderen Methoden oft nur sehr umständlich 
zu lösen sind. Um die Laplace-Transformation aber 
schnell und sicher anwenden zu können, braucht man 
ein möglichst umfangreiches Verzeichnis von Funktio- 
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nen (Oberfunktionen) mit den dazu gehörenden La- 
place-Transformierten (Unterfunktionen). ee 
das vorliegende Buch, das die klassische e; 

nale Transformation behandelt, in einem Umfang, wie 
er bisher in Deutschland noch nicht gegeben wurde. 
Im zweiten Teil des Buches finden sich fast 800 solche 
Korrespondenzen, die nach dem Gesichtspunkt der 
praktischen Brauchbarkeit ausgewählt sind. Die An- 
ordnung ist nicht die bisher übliche nach Oberfunktio- 


nen, sondern Verfasser ordnet nach Unterfunktionen, _ 


Das ist deswegen vorteilhaft, weil man bei der Lösung 
von Differentialgleichungen zu einer Transformierten 
die Ausgangsfunktion zu finden hat. Übrigens ist am 
Schluß dieses zweiten Teiles ein Register der auftreten- 
den Oberfunktionen gegeben, das die Funktionssym- 
bole, Definition und an der Funktionen und Ver- 
weise auf die Abschnitte und Nummern, in denen sie 
als Oberfunktionen vorkommen, enthält. 

Um nun aber das Buch auch für solche Benutzer 
brauchbar zu machen, die sich nicht eingehender mit 
der Theorie der Laplace-Transformation beschäftigt 
haben, schickt der Verfasser den Tabellen einen ersten 
Teil mit einer Anleitung zu ihrem Gebrauch voraus. 
Dieseist außerordentlich lesbar geschrieben. Sie bringt 
nach Definition der Transformation und Angabe der 
benutzten Bezeichnungen in übersichtlicher Zusam- 
menstellung die Haupteigenschaften und die wichtig- 
sten Regeln zum Rechnen mit diesen Transformationen. 
Jeder Leser kann so unter Benutzung dieser Regeln die 
Transformation praktisch anwenden, vor allem, wenn 
er die Beispiele Aurhgensheiie hat, an denen ihr Ge- 
brauch zur Lösung gewöhnlicher und partieller Diffe- 
rentialgleichungen erläutert wird. In diesen Beispielen 
wird auch ein Fall vorgeführt, wo die gedankenlose Ver- 
wendung der Umkehrformel zu falschen Resultaten 


NACHRICHTEN E 


Max Planck }. 


Am 4.10. 1947 ist Max Planck, der Schöpfer 
des Wirkungsquantums, in Göttingen an einer Lungen- 
entzündung gestorben. Es gibt kein Land der Erde, 
dessen Physiker nicht des Mannes gedächten, der durch 
sein Strahlungsgesetz eine neue Epoche in der Physik 
eingeleitet hat, und die ihn nicht in der Geschichte der 
Physik an die Seite vonNewton, Maxwellu.a. 
stellen. Viele Arbeiten aus dem Gebiete-der Thermo- 
dynamik, die ihn schon in jungen Jahren beschäftigte, 
und die sich mit den Fragen der Entropie befaßten, 
sind von ihm verfaßt worden. Die grundlegendste, die 
das Weltbild der Physik zu verändern vermochte, ist 
im Jahre 1900 während seiner Lehrtätigkeit an der 
Berliner Universität, an der er bis zu seiner Emeritie- 
rung wirkte, entstanden. Sie ist im Anschluß an Dis- 
kussionen über das Wiensche Strahlungsgesetz in 
einer Sitzung der Deutschen Physikalischen Gesell- 
schaft vorgetragen worden. Versuche, klassisch das 
Wiensche Strahlungsgesetz zu erklären, führten: zu 
eklatantem Mißerfolg. Plan-ck hat dann den Ge- 
danken ausgesprochen, daß sich diese Schwierigkeiten 
durch Annahme diskreter Energiequanten beseitigen 
lasse. Er konnte so die Existenz des Wirkungsquan- 
tums begründen. Diese Plancksche Entdeckung 
eröffnete eine neue Epoche in der Physik und hat vor 
allem die Physik der atomaren Vorgänge befruchtet. 
Es gibt darin kaum einen Zweig, in dem nicht das Wir- 
kungsquantum eine entscheidende Rolle spielt. 

Neben seiner Leistung war es vor allem sein vor- 
nehmer Charakter, der ihn befähigte, noch in seinen 


letzten Jahren die Geltung der Wissenschaft in dem 
Deutschland der Hitlerzeit zu verteidigen. 


München: Der Privatdozent der Universität Tü- 
bingen, Dr. Fritz Bopp, wurde zum a. o. Prof. für 
theoretische Physik und angewandte Mathematik an 
der Universität München ernannt. 


Bonn: Der apl.a. o. Prof. der Universität Heidelberg, 
Dr. Hilmar Wendt. liestin diesem Winter als Gast- 
professor an der geodätisch-kulturtechnischen Ab- 
teilung der landwirtschaftlichen Fakultät der 
Universität Bonn. 


Frankfurt .a.M.: Am 30. 6. 1947 wurde Herrn 

Dr. G. Lorentz eine Dozentur für Mathematik an 
der Universität Frankfurt übertragen. 
Am 10. 10. 1947 wurde Fräulein Dr. R.Moufang 
mit der kommissarischen Wahrnehmung des Extra- 
ordinariates für Mathematik an der Universität 
Frankfurt beauftragt und am 19. 12, 1947 zum 
apl. Professor ernannt, j 


Berichtigung: 


Die an letzter Stelle unter den Nachrichten des Hef- 
tes 5/6 gebrachte Notiz bedarf der Berichtigung. Es 
muß heißen: 

Zürich: Prof, Stephen Timoshenko an der 
Stanford University Cal. USA, wurde zum Dr. sc. 
techn. h. c. der Eidgenössischen Technischen Hoch- 
schule Zürich gewählt. 
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